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Zusammenfassung 



Das Axion ist ein pseudo-Nambu-Goldstone Boson. Es tritt in Erscheinung nach 
dem spontanen Bruch der Peccei-Quinn Symmetrie, die als Losung des starken 
CP-Problems vorgeschlagen wurde. Andere pseudo-Nambu-Goldstone Bosonen, 
postuliert in einigen Erweiterungen des Standardmodells, werden "Axion-Like Par- 
ticles" (ALPs) genannt, wenn sie bestimmte Eigenschaften mit dem Axion teilen, 
insbesondere die Kopplung an zwei Photonen. Bis jetzt waren alle Suchen nach A- 
xionen und ALPs erfolglos. Dies bedeutet, dass deren Kopplungen extrem schwach 
sein miissen. Allerdings konnen Axionen und ALPs einige beobachtbare astro- 
physikalische und kosmologische Auswirkungen haben, anhand derer man den Pa- 
rameterraum dieser Teilchen einschranken kann. 

Wir konzentrieren uns auf Schranken aus der Kosmologie, die ein ideales Feld 
fur die Untersuchung von Axionen und ALPs darstellt. Insbesondere untersuchen 
wir als erstes die Moglichkeit einer Axion- und ALP-Population, die wahrend 
der fruhesten Augenblicke des Universums entstanden ist. Die Bedeutung dieser 
Analyse ruhrt daher, dass Axionen und ALPs wegen ihrer schwachen Wechsel- 
wirkung und der besonderen Produktionsmechanismen ideale Kandidaten fur die 
dunkle Materie sind. Schliefilich betrachten wir die Folgen des Zerfalls dieser 
Teilchen fiir bestimmte kosmologische Observablen, namlich fiir das Photonen- 
spektrum von Galaxien, fiir den kosmoschen Mikrowellenhintergrund, fiir die ef- 
fektive Zahl an Neutrinos und die urspriingliche Haufigkeit der Elemente. Unsere 
Schranken stellen den striktesten Test eines fruhen Zerfalls von Axionen und ALPs 
dar und schliefien einen Teil des Parameterraums von ALPs aus, der ansonsten ex- 
perimentell nur schwer zuganglich ist. 



Abstract 



The axion is a pseudo-Nambu-Goldstone boson. It appears after the spontaneous 
breaking of the Peccei-Quinn symmetry, which was proposed to solve the strong- 
CP problem. Other pseudo-Nambu-Goldstone bosons, postulated in some ex- 
tensions of the standard model of particle physics, are called axion-like particles 
(ALPs) if they share certain characteristics with the axion, in particular a coupling 
to two photons. Thus far, axion and ALP searches have been unsuccessful, indi- 
cating that their couplings have to be extremely weak. However, axions and ALPs 
could be responsible for some observable effects in astrophysics and cosmology, 
which can also be exploited to constrain the parameter space of these particles. 

We focus on limits coming from cosmology, which is an optimal field for study- 
ing axions and ALPs. In particular, we first investigate the possibility of a pri- 
mordial population of axions and ALPs arising during the earliest epochs of the 
universe. The importance of this analysis lies on the fact that axions and ALPs are 
ideal dark matter candidates because of their faint interactions and their peculiar 
production mechanisms. Finally, we consider the consequences of the decay of such 
a population on specific cosmological observables, namely the photon spectrum 
of galaxies, the cosmic microwave background, the effective number of neutrino 
species, and the abundance of primordial elements. Our bounds constitute the 
most stringent probes of early decays and exclude a part of the ALP parameter 
space that is otherwise very difficult to test experimentally. 
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Preface 



Day after day, cosmology has become more useful as a tool for particle physics. In 
its earliest epoch, the universe provided the environment with the highest known 
energy, much higher than those presently tested in colliders. We still miss a particle 
physics description of such high energy scales, therefore many open cosmological 
questions are waiting for a satisfying answer in terms of microphysics. Every 
proposal of physics beyond the standard model of particle physics (SM) must 
deal with cosmology, notably providing mechanisms for inflation and baryogenesis 
and candidates for dark matter and dark energy. All these issues provide further 
motivations to reach a complete picture of particle theory. 

On the other hand, through the study of those phases of cosmic evolution, when 
both the microphysics and cosmological descriptions are well settled, physicists 
can put severe constraints on physics beyond the SM. Big-bang nucleosynthesis 
(BBN), the cosmic microwave background (CMB) and large scale structure (LSS) 
give precise indications about what happened when the universe was older than a 
second. Moreover, in the last few years the vastly increasing amount and quality 
of data has propelled cosmology to its precision era. New particle physics models 
have to face these broad data sets and must be consistent with them. Cosmology 
is therefore a crucial testing ground for particle physics. 

But cosmology does not play only the killjoy role. In particular, because of the 
intrinsic long exposure time which is peculiar of cosmological phenomena, together 
with the high luminosity provided by the universe itself, cosmology gives the pos- 
sibility of testing very weakly interacting particles otherwise not accessible to high 
energy particle physics experiments. Axion and axion-like particles (ALPs) are 
among these species, and they are the main characters of this dissertation. Very 
stringent limits on their mass and coupling constants come from observations of 
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high-intensity astrophysical settings like globular clusters or the sun. Cosmology 
also plays an important role in the quest for finding them, and, as we will see, the 
bounds that it provides are complementary to those provided by astrophysics. 

In the first chapter, we introduce the QCD axion and ALPs, providing some 
theoretical motivation. We also review briefly the present status of the exploration 
of their parameter space through non-cosmological approaches. We already men- 
tioned the astrophysical observations, and here we describe the main experimental 
approaches to axion and ALP physics. We conclude this chapter with some hints 
for the presence of such particles, supported by astrophysical clues. 

The primordial production of a population of axions and ALPs is the topic of 
chapter 2. This is the basic question to scrutinise before even to start the discussion 
about the bounds cosmology could provide: without a relic population there can 
be no cosmological limits. Moreover, establishing a relic population of axions or 
ALPs is an intriguing topic in its own right, as it could provide the solution to the 
dark matter mystery. The analysis of the stability of these particles is therefore 
the optimal way to conclude this chapter, and to introduce the core question of 
this dissertation. 

If axions and ALPs are unstable, they decay into photons. Their decay products 
could influence the evolution of the universe and the cosmological observables. In 
chapter 3 we treat the limits coming from the late decays of axions and ALPs. 
These are obtained studying how the decay affects the CMB and considering the 
possibility of directly detecting the photons emitted in the spectrum of galaxies 
and the extragalactic photon background. 

At earlier times, when the universe was hot and dense enough to rapidly lose the 
direct imprint of the decay products via thermal scattering, the decay would have 
a more subtle and indirect influence. In chapter 4 we notice that the decay, which 
injects a large amount of photons in the primordial plasma, can effectively dilute 
baryons and neutrinos, if it happens after about ten milliseconds since the big 
bang. The primordial elemental abundances and the amount of radiation during 
the first minutes of our universe are strongly affected by this dilution and provide 
further limits on the axion and ALP parameters. To conclude this chapter, we 
reanalyse the constraints we found considering other decay channels besides the 
photon one. Finally, chapter 5 is devoted to summarise our arguments and to draw 
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our conclusions. 

The research on axions and ALPs will be one of the main frontiers of particle 
physics in the near future, since these pseudoscalars can solve some of the unre- 
solved problems of particle physics. Their experimental discovery would be a true 
milestone along our path to understand Nature. The technological challenge to 
reach this aim is pushing our capability of measuring extraordinarily small signals, 
but a long path remains to be covered. Astrophysical and cosmological observa- 
tions provide directions about the particle parameters, offering a guidance in the 
design of devoted discovery experiments. The improving quality and quantity of 
astrophysical and cosmological data gives us confidence on the possibilities the sky 
offers. That is why we consider this kind of analysis of fundamental importance 
in the quest for new physics. 



CHAPTER 1 



Axions and axion-like particles 



1.1 The strong-CP problem and the axion 

Quantum chromodynamics (QCD), the theory of strong interactions, includes a P, 
T and thus CP violating term [4-7], 

C* = 0^^G?G>?, (1.1) 

where G? a is the gluon field strength, a s the fine structure constant for colour 
interactions and the so called 9 angle, — n < 9 < ir, is the effective parameter 
controlling CP violation. Here, and in the following, we will use natural units with 
H — c — kj$ — 1, where k B is the Boltzmann constant. It is common practice to 
write the dual tensor e M j ypo -G^ <T /2 simply as G afiv , thus we will simply write GG for 
the trace in the Lagrangian (1.1). The sources of this 0-term are the Adler-Bell- 
Jackiw anomaly of the axial current in QCD [8, 9] and the topology of the QCD 
vacuum [10]. There are no theoretical hints about which value between — n and 
+7T the parameter 9 could choose, so we just expect it to be an 0(1) quantity. 

The gluon field transforms like a 4- vector under C, P and T. We can define 
the coloured electric and magnetic fields as 

E k = G ok = d o A k _ d k A o + g j abcA o A k s (L2a) 

B k a = e ijk G ija = e ijk (d t A ja - & A ia + g s f abc A ib A jc ) , (1.2b) 

and they transform like their electromagnetic counterparts under C, P and T. In 
the definitions (1.2), g s is the strong interaction coupling constant and f a b c are 
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the SU(3) structure constants. The colour indices are abc, while ijk are spatial 
indices. Because of the properties of the completely antisymmetric tensor e^ pcr , 
the Lagrangian (1.1) can be written as GG = —AB a ■ E a . The scalar product of a 
polar vector and an axial vector violates P, T and thus CP, once CPT is taken for 
granted. 

Whether CP is a good symmetry for QCD is not a fundamental quastion but a 
phenomenological one. At low energy the Lagrangian (1.1) induces electric dipole 
moments in baryons, which have not yet been observed. In particular, many 
measurements have been performed on the neutron one. Calculations predict the 
neutron electric dipole moment to be d n ~ 10" 16 ^ecm [11,12], where e is the 
electron electric charge. Currently, the best experimental limit is d n < 0.29 x 
10~ 25 ecm [13], which translates into 

\e\ < i(r 10 . (i.3) 

The essentials of the strong-CP problem are all here: we were not expecting such 
an extremely small value for 9. 

Roberto Peccei and Helen Quinn proposed an elegant solution to this puz- 
zle, introducing a global chiral U(l) symmetry, called the Peccei-Quinn symmetry 
£7(1)pq [14,15]. This symmetry is spontaneously broken, and its Nambu-Goldstone 
boson (NGB) a, the celebrated axion [16,17], couples to gluons because the £/(1) P q 
symmetry is violated by the colour anomaly. Thus, another term involving GG 
enters in the QCD Lagrangian, 

where f a is the order parameter associated with the breaking of [7(1) P q. For the 
moment we ignore the subtleties related to the definition of f a . At energies below 
the confinement scale of colour interactions, Aq CD ~ 1 GeV, gluons and quarks 
have to be integrated out, and QCD is described by an effective chiral Lagrangian. 
The terms (1.1) and (1.4), together with the anomalous contribution of the U(l) A 
symmetry of the effective chiral QCD Lagrangian, become an effective potential 
that can be parametrised as 



V(a) ~ A 4 QCD 



cos I — - — - 9 

J a Jrj 



;i.5) 
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Here, rj' is the pseudo-NGB (PNGB) of the anomalous U (1)a, and f v is a parameter 
with dimensions of a mass. The periodicity of the potential is due to the topological 
and instantonic nature of the GG Lagrangian. The configuration of minimum 
energy for the potential (1.5) is realised by the CP-conserving linear combination 
(a) — (9 — (rf')/f w ) /„, as required by the Vafa-Witten theorem [18]. Since the a 
field is massless, it can align itself with a null energy cost to the CP-conserving 
point: the strong- CP problem is solved, independently of the value of f a , 9 and (77'). 

Because of the effective potential (1.5), the mass of the eigenstate rj' + af v /f a 
is rrirji ~ Aq CD // t; , which solves Weinberg's U(l) problem 1 . The orthogonal com- 
bination, which mainly consists of axion, mixes with the other scalar mesons in 
the particle spectrum and acquires mass [20]. It is this combination which is usu- 
ally called axion, and whose phenomenology is studied. In first approximation - 
considering only up and down quark contributions — the axion mass is [21] 

m w f w y/m~m2 /10 6 GeV\ 
m a = — ■ ~ 6 eV , 1.6 

f a m u + m d V fa J 

where m n = 135 MeV is the neutral pion mass and f n = 92 MeV is the pion 
decay constant. The errors in the measurements of the light quark mass ratio 
establish the range m u /m a ' = 0.3-0.6 for this quantity [22], which translates into 
a 10% uncertainty on the axion mass, m a = 5.23-6.01 eV for f a = 10 6 GeV. The 
preferred value m u /md = 0.56 gives us the result (1.6) after rounding it to 10% 
accuracy. 

In the original Peccei-Quinn model, the additional chiral symmetry was im- 
posed on the SM through two Higgs doublets, linking the spontaneous breaking 
of L r (l)pQ to the electro-weak symmetry breaking, thus f a ~ -Eew — 246 GeV. 
However, this first attempt to solve the strong-CP problem through the Peccei- 
Quinn symmetry was quickly ruled out, because the axion did not show up in 
experimental data [23]. 



1 An approximate U(l) axial- vector current would require the presence of a pseudoscalar boson 
with mass smaller than ^/Sm^, which is not observed [19]. This fact was considered a problem 
before realising that the axial current is violated by colour anomaly [8,9]. Thus the #-term in 
the QCD Lagrangian solves the U(l) problem, and this provides significance to the 9- vacuum. 
Of course, the presence of the axion is not required to solve the U(l) problem, but it helps to 
recover the observed CP conservation. 
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A way to save the Peccei-Quinn mechanism, bypassing the experimental limits, 
is to raise the parameter f a by several orders of magnitude. The scale f a suppresses 
both the axion mass and couplings, a fact that will become more clear shortly. For 
the moment it is sufficient to notice that since the axion is usually introduced as a 
phase in the Higgs sector, it needs to be normalised by a constant with dimension 
of energy. It is therefore the dimensionless ratio (a// ) that always appears in 
effective Lagrangians for axions. Increasing f a automatically lowers the axion 
couplings to SM particles, leaving untouched the validity of the solution to the 
strong-CP problem, as we have seen before. There are several implementations 
of this idea, which are called invisible axion models. The first of them was the 
Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion model [24,25], which we discuss 
here in more detail for its simplicity. A Dirac quark field Q — a colour triplet 
in the fundamental representation with no bare mass — and a SM gauge group 
singlet Higgs-scalar S are added to the standard model. The Lagrangian we need 
is 



where 7 M are the Dirac matrices, g s is the strong interaction coupling constant, 
the gluon field, A a are the Gell-Mann matrices, Qur) is the left (right) handed 
projection of the Dirac spinor, y the Yukawa coupling, A > and v a 3> -Eew 
are the parameters of the Mexican hat potential, the last one with dimensions of 
energy. For the moment we do not take into account the interaction of Q with the 
electroweak gauge fields. Under the global chiral C/(1)pq, the fields transform as 



£ksvz 



l -a 
- 7 GT + O-r-GG + tQYd.Q + g s A^ a Q lfl X a Q 

4 07T 




(1.7) 



Qr ->■ Que 



,ia/2 



(1.8a) 
(1.8b) 
(1.8c) 



-ia/2 



S ^ Se 



leaving the Lagrangian £rsvz invariant at the classical level. 

It is convenient to express the scalar Higgs singlet in its polar form, 




(1.9) 
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Figure 1.1: Triangle loop diagram for the effective axion-gluon interaction of 
equation (1.12). 



At energies lower than v a , S rolls toward the minimum of the Mexican hat po- 
tential, and we can make the substitutions p — > v a . If we limit our physics consi- 
derations to energies lower than the singlet mass, we can keep only the phase field 
a. This field changes as a — > a + av a under PQ transformations (1.8). However, 
because of the invariance of the Lagrangian (1.7), the energy cost of this shift 
is null, and therefore a is the NGB for {7(1)pq. Under this approximation, the 
Lagrangian is now 

1 -a 
£'ksvz = - ~ A G^G a , u + 9-^GG + iQfd^Q + g s G^ a Q lfl \ a Q 

- yv a (qU^Qr + Q R e'^Q L ) + . (1.10) 

The spontaneous breaking of the Peccei-Quinn symmetry produces an effective 
mass term uiq = yv a for the quark field, once the phase terms between the left 
and right spinors are removed. To achieve this, we chiral rotate the quark field, 

This transformation adds to the Lagrangian (1.10) the terms 

5C = -^Q l5 YQ-^-GG, (1.12) 

2v a 87T v a 

where the last term is the contribution of the colour anomaly, originating from the 
triangle diagram of figure 1.1. Dealing with energy scales much lower than itlq, 
the quark field can be integrated out, and the remaining terms of the Lagrangian 
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are 

£ksvz 



-\G^G ai ^ + ($ - ^GG + l -d,ad»a , (1.13) 

among which there is the C aG Q term of equation (1.4). 

If many heavy quark fields Q are introduced in this model, each of them 
would produce a contribution like equation (1.12). Under a PQ-transformation, 
these heavy quarks transform if a PQ charge qpg is assigned to them, i.e. Q 3 L — > 
Q J L exp(iqpQa/2) and Q R — > Q R exp(—iq P Qa/2) if S — > Se ta . This leaves the La- 
grangian invariant under the condition gpq = 1. In the effective Lagrangian we 
define the Peccei-Quinn scale or axion decay constant to be 

fa = ]±, (1-14) 

where Af = Ylj Qpq counts the number of species that are PQ-charged, and the 
form of C aG Q of equation (1.4) is restored. 

A triangle diagram, like that of figure 1.1, can create other anomalous contri- 
butions to equation (1.12), if the heavy quark fields are also coupled to electroweak 
gauge bosons. In the low-energy regime, this means the presence of a two-photon 
coupling for the axion, 

^77 = -^aFF , (1.15) 

f" being the photon field strength and F^ v its dual. Anyway, because of axion- 
meson mixing, the two-photon coupling arises even if the heavy quark field is 
completely decoupled from the SU{2) x U(l) gauge sector. The coupling constant 
g ail is actually the sum of two contributions [5], 



a 
2vr/ ( 



S£ (^E«pq (ta.) - - 3l ^: j = ^j-c: . (Me) 



where a is the fine structure constant. The first term is the electromagnetic anoma- 
lous contribution of the heavy quarks, q P q ^ EM ^ being the Peccei-Quinn (electro- 
magnetic) charge of the i-th Q field, and the second one is produced by the axion- 
meson mixing. This model dependence can be parametrised with the dimensionless 
quantity C 7 . In the literature, the KSVZ axion is usually defined to have no elec- 
tromagnetic anomaly and thus to have |C 7 | ~ 1.95 for m u /md = 0.56. Because of 
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the uncertain up and down quark mass ratio, the parameter |C 7 | could be anywhere 
between 1.92 and 2.20. 

Another popular axion model was proposed by Zhitnitsky [26] and by Dine, 
Fischler and Srednicki [27] and is called the DFSZ axion. Again, the axion is 
related to the phase of a Higgs-like scalar singlet S which is added to a two Higgs 
doublet extension of the SM. The scalar S is coupled only to the Higgs doublets 
H u and H d in the Higgs potential, 

V(H U , H d , S) = \ u (HjH u -v 2 u ) 2 + \ d (HjH d - v 2 ) 2 + A (S*S - v 2 a ) 2 
+ X {HjH u ) (HjH d ) + C {HjH d ) (HjH u ) 

+ [ lu (HjH u ) + ld (HjH d )] S*S + 7 [(HjH d ) S 2 + ha] , (1.17) 

where X u , X d , A, \i C> lw Id and 7 are real dimensionless parameters of the 
potential. This potential is invariant under the PQ-transformation 

H u ^H u e iX \ H d ^H d e iXd , S ->■ Se i{x ^- Xd)/2 . (1.18) 

The model can equally work if the last term of the potential is 7 ^(H^Hd) S + h.c] , 
but this time 7 has to have mass dimensions and the PQ-transformations have to 
be adapted. The Yukawa Lagrangian is 

C Y = yue a bQ a L Hl b u R + y d Q L H d d R + y d L L H d l R + h.c. , (1.19) 

where e is the 2x2 antisymmetric matrix and, in this case, Ql and Ll are 
the left-handed SU(2) quark and lepton doublets of the SM. The right-handed 
charged lepton component is l R , and the Yukawa couplings are the ys. The vacuum 
expectation values (VEV) of the Higgs fields can be written as 

The two Higgs doublets have four degrees of freedom each, and three of them 
are SU(2) phases, while S has got just a radial excitation and a phase, like in 
the KSVZ model. For phenomenological reasons we want v a >■ v — ^Jv\ + v\ ~ 
E EW . Once C/(1) P q is broken, the radial part of S freezes on its VEV, while its 
phase degree of freedom remains free. At energies below the EW phase transition, 
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the two Higgs doublets can be written in the unitary gauge deleting from the 
theory the redundant massless degrees of freedom, which are "eaten" by the gauge 
bosons. At this point, the neutral orthogonal combination of degrees of freedom 
which survives mixes with the phase of the S boson due to the last term of the 
potential (1.17). In this mixing, one combination is massive and is the Higgs boson, 
while the orthogonal massless NGB is the axion a. Redefining the Higgs doublets 
through transformations like (1.18) to reabsorb the axion phase in the last term 
of the potential (1.17), makes a reappear in the Yukawa terms (1.20). Again, 
the axion phase can be reabsorbed with some chiral rotations of the fermionic 
fields, producing a set of anomalous terms like equations (1.12) and (1.15), but 
this time with the SM fermions instead of the heavy quark Q. In the DFSZ model, 
the derivative couplings in (1.12) play a role in the low energy theory, since they 
involve the light fermions too. The strengths of the couplings depend on the PQ 
charges X u and X d , and are therefore model dependent. 

We presently do not know which are the features of the high energy model 
which gives origin to the axion. However, we can adhere to the proposal of Georgi, 
Kaplan and Randall [28], and write a generic low energy effective theory. In the 
range between E EW and Aq CD , the Lagrangian is 

C a = \d»ad»a -j^GG- 9 -^aFF - ?g £ C^ /757 <% . (1.21) 

All the model dependencies are hidden in the coupling coefficients g ail and Cf, 
where the index / runs over all the SM fermions. The definition of the axion- 
photon coupling constant g ary is provided by equation (1.16). In particular, the 
KSVZ model predicts Cf — for all the leptons and the ordinary quarks at tree 
level. In the DFSZ model, the coupling coefficient to electrons is C e = cos 2 (j3) /3, 
while the couplings to the up and down quarks are respectively C u = sin 2 (/3) /3 
and Cd = cos 2 (/3) /3, where /3 is the ratio of the vacuum expectation values of the 
two H u and fields [22]. Below Aqcd, gluons and quarks confine so we have to 
write an effective Lagrangian including the couplings to nucleons and mesons, 



A. = ^ad^a - m 2 Jl 



cog j _ g^ aFp _ |a ^ C/ ^- /75yV/ + Cc 



(1.22) 



1.2 Enlarging the parameter space: axion-like particles 
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The potential is the same of equation (1.5) up to a reparametrisation. This time, 
the sum over / covers the light SM leptons and the nucleons. The coupling co- 
efficients for axion-proton and axion-neutron interactions are C v = —0.55 and 
C n = 0.14 for m u /rrid = 0.3 and C p = —0.37 and C n = —0.05 for m u jm^ = 0.6 
in the KSVZ model [22]. The DFSZ axion has couplings to nucleons of the same 
order which depend also on f3 [22] . Among the interactions with mesons, only the 
pion-axion one is interesting for our purposes. Therefore in Lagrangran (1.22) we 
express C a7T as [22] 

C a7T = ^f-^JL (f + 7T°7r-^7r + - 27r + 7T-^7T ) , (1.23) 

J a J it 

where C n is again a model dependent constant. 



1.2 Enlarging the parameter space: 
axion-like particles 

As explained in the previous section, several axion models exist, each of them sol- 
ving the strong CP problem, but providing different couplings to the SM particles. 
Moreover, the PQ-mechanism works for every value of f a , and the clues about 
the characteristic PQ-scale are only speculative. The experimental search for the 
axion requires therefore not only to explore many orders of magnitude in f a or 
m a , but even to scan the model dependencies held in the coupling coefficients of 
Lagrangians (1.21) and (1.22). 

After the seminal paper by Sikivie [29], the most relevant axion direct searches 
try to exploit the two-photon coupling. However, as we will show in the next 
section, the discovery task is tough, the axion being very weekly coupled for the 
allowed range of f a . Anyway, this kind of experiments could find in principle any 
kind of pseudoscalar particle <fi coupled to photons through a term like 

C = -Ig^F^ , (1.24) 

which mimics the axion-photon interaction (1.15) [30,31]. Therefore, the interest 
on such hypothetical particles that could couple to two photons has risen in the 
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last years, and the name axion-like particle (ALP) has been coined for them. These 
ALPs can in principle also couple to other particles besides photons, however we 
will always deal with the interaction Lagrangian (1.24), as it is typically the most 
important one for the phenomenology at low energies. We will comment about 
other possible couplings at the end of chapter 4. 

ALPs are even more interesting on the theoretical side, as they can arise in the 
low energy spectrum of many extensions of the SM. An ALP can appear in a theory 
as a PNGB of a continuous global symmetry. Examples of these symmetries are 
related to particle flavour [32], lepton number [33,34] or the i?-symmetry in super- 
symmetry [35,36]. When a continuous global symmetry is spontaneously broken, 
massless NGBs appear in the low energy theory as phases of the high energy degrees 
of freedom. Since phases are dimensionless, the canonically normalised theory at 
low energies always involves the combination 0//^, where is the NGB field and 

is a scale close to the spontaneous symmetry breaking (SSB) scale. The ALP 
could be for instance related to the generation of right-handed neutrino masses, 
and hence have a decay constant at an intermediate scale, like ~ 10 10 — 10 12 GeV; 
alternatively it could be associated with a grand-unification theory (GUT), and 
have a decay constant at the corresponding scale f<p ~ 10 15 GeV. From the no- 
hair theorem, we know up to some extent that black-hole dynamics violates global 
symmetry conservation. Therefore unbroken global symmetries can not exist in 
theories with gravity, and we should have PNGBs instead of NGBs. There are 
many possibilities for breaking the shift symmetry besides gravity effects, explici- 
tly or spontaneously, perturbatively or non-perturbatively. 

Moreover, the observation that in string theory ALPs appear in all compactifi- 
cations has raised even more attention to them [37]. These so-called string axions 
share the NGB properties (having a shift symmetry and being periodic) but with 
the natural size of being the string scale. 

The fact that the ALP acquires a mass implies that in the model Lagrangian 
a potential has to be included. Taking inspiration from the axion case, the ALP 
potential can typically be parametrized as 



If the dynamics explicitly breaking the associated global continuous symmetry have 




(1.25) 
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a characteristic scale A, the ALP mass is parametrically small, since it is suppressed 
by powers of A/f^. The phenomenology of the SM requires A to be related to 
physics beyond the electroweak scale, i.e. A > TeV (which implies ^> m a 
for = f a ), or to belong to a hidden sector. We have here no preconceptions 
regarding ALP mass, for it depends on the unknown ratio of two unconstrained 
energy scales. 

The dimensionful coupling parameter in equation (1.24) can be parametrised 

as 

<L26) 

In the simplest is an integer, but this is not true in general when the 

ALP mixes, either kinetically or via symmetry-breaking effects with other ALPs 
or with pseudoscalar mesons. For string axions, the coupling to photons is related 
via a loop factor to either the string scale or the Planck scale, or it could be even 
weaker, so an ALP with a large coupling would restrict the string scale to be low. 
For string and field theoretical models the most interesting values are therefore 
gj, ~ Kr n -l(r 15 GeV-\ ~ KT 19 GeV" 1 and ~ KT 21 GeV" 1 corresponding to 
intermediate, GUT or Planck scales. 

We will adhere to the phenomenological approach to leave the ALP parameters 
rriff, and <70 free to span many orders of magnitude, exploring the consequences of 
the presence of these particles in order to limit the ALP parameter space, or to 
find hints of their existence. 



1.3 A compendium of limits 

The search for axions and ALPs has not yet been successful. Up to now we have 
only indications about where these scalars are not and some hints about where 
they could hide. The present bounds on the axion mass, and consequently on the 
PQ-scale, are plotted in figure 1.2. These bounds are indicative, as the couplings 
can change according to the axion model, but they give a picture of the situation, 
especially considering 0(1) coupling coefficients. The red and brown bounds in 
the first line, labelled Cold DM, Topological defect decay and Hot DM are 
the bounds provided by cosmology and we will thoroughly discuss them in the next 
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Figure 1.2: Axion limits. The red bounds in the first line come from cosmological 
considerations. The blue ones in the second line are due to direct measurements. 
The green limits refer to astrophysical arguments related to stellar evolution. 
These bounds are described in more detail in the text. The dotted ellipse suggest 
the values of m a that could help fitting recent white dwarf cooling data. 



chapters. The blue bounds in the second line come from direct measurements of 
astrophysical and cosmological quantities and from experiments, while the green 
ones are related to astrophysics. 

As we have already mentioned, very soon after it was proposed the axion was 
ruled out in the keV-MeV mass range. Measurements of heavy quarkonium state 
decays or of nuclear de-excitations and beam dump or reactor experiments have 
found no evidence for an axion coupled to fermions or nucleons up to f a > 10 4 GeV, 
which in term of mass is m a < 0.6 keV [7]. Pure particle physics experiments 
can only test values f a in a rather low energy range. The bound they put dies 
out when f a is high enough to fade the production and detection probabilities 
through the suppression of the axion couplings. In figure 1.2 this bound is labelled 
Laboratory. 

The Telescope region, m a = 3-27 eV and f a = 2.3 x 10 5 -2.1 x 10 6 GeV, is 
excluded by the non-observation of photons that could be related to the relic ax- 
ion decay a — > 77 in the spectrum of galaxies and in the extragalactic background 
light [38-41]. The Axion Dark Matter eXperiment (ADMX) has provided some 
constraints in the region m a = 1.9-3.5 /ieV, which translates into f a = 1.8 x 10 12 - 
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3.3 x 10 12 GeV, after the sensitivity to test the axion was almost reached [42]. This 
collaboration is looking for axion dark matter (DM) using an haloscope. We will 
describe this instrument in the next chapter. The plans of the ADMX collabora- 
tion are to scan the mass range enclosed in the dashed region called ADMX in 
figure 1.2 after some upgrades [43-45]. The two latter limits are closely related to 
cosmology and to the axion being a DM component, we will therefore refer again 
to them in the following chapters. 

Stars represent a prolific environment for the production of light and weakly- 
coupled particles, as we have learned for example in the neutrino case [46]. The sun 
is surely the brightest axion source in the sky. A photon can convert into an axion 
if it interacts with an external magnetic or electric field by means of the two-photon 
coupling of equation (1.15) as in figure 1.3 [47]. This is the so called Primakoff 
effect, which was first proposed for the creation of mesons in the electric field of 
nuclei [48]. Therefore, the high concentration of thermal photons, together with 
the strong electromagnetic fields of the stellar plasma, makes the sun, and stars 
in general, a rich soil for axion and ALP production. The total axion luminosity, 
calculated using the standard solar model, is L a ~ (g ary /10~ 10 GeV _1 ) 2 lO~ 3 L , 
where L & = 3.90 x 10 25 W is the solar luminosity in photons [49]. To detect the flux 
of axions, several solar axion telescopes, like SUMICO [50] and CAST [51], have 
been built. These helioscopes are essentially vacuum pipes. They are permeated 
by a strong magnetic field to exploit the inverse Primakoff effect to convert axions 
back into photons [29]. Axions can enter the telescope because of their very weak 
interaction with matter, and successively be detected once they have oscillated into 
photons. Among the helioscopes, CAST currently gives the strongest constraints: 
its results exclude the part of the axion parameter space, m a = 0.39-0.64 eV and 
f a = 9.8 x 10 6 1.6 x 10 7 GeV, labelled CAST. Recently, a new proposal for an 
axion helioscope has appeared, the International AXion Observatory (IAXO) [52]. 
The hope is to improve the sensitivity to g ail of at least one order of magnitude 
with respect to CAST and therefore, in the most optimistic scenario, to explore 
the area labelled IAXO which is enclosed by the blue dashed line. 

Waiting for IAXO and its results, the best upper limits on the axion mass in 
figure 1.2 come from stellar evolution arguments. A globular cluster is a gravi- 
tationally bound system of stars. The stars belonging to a globular cluster have 
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Figure 1.3: Diagram of the Primakoff effect. 

approximately the same age and they differ only in their initial mass. Since the 
more massive a star is, the faster it evolves, a globular cluster gives the possibility 
to study a broad sample of stellar evolution stages and to estimate how long each 
phase lasts. In particular, if axions are produced inside a star and escape, they 
provide an additional cooling channel, besides the photon and neutrino ones. If 
there are more efficient energy release channels, the nuclear fuel consumption has 
to be faster, and thus the ageing quicker. Counting the stars in each evolution 
stage inside a globular cluster permits us to study how fast the fuel consumption 
is and therefore to put bounds on the production of axions in stellar cores. The 
best constraints come from the stars which have reached the helium burning phase, 
which are called horizontal branch (HB) stars because of the position they occupy 
in the Hertzsprung-Russel diagram. The non-standard energy loss prolongs the 
red giant (RG) phase and shortens the HB one [46]. Counting the RG and the HB 
stars in globular clusters and comparing the two numbers it is possible to evaluate 
the axion production rate in stars, and to obtain the two HB Stars bounds in 
figure 1.2. In particular, if the axion is directly coupled to the electron, i.e. C e is 
(9(1), it has a significant production channel more which is reflected in the broader 
exclusion bound. HB stars have a typical core temperature of T ~ 10 8 K ~ 10 keV. 
The thermal distribution of photons, averaged over the large volume of the star, 
still includes many 7s that are energetic enough to efficiently produce axions if 
their mass is not m a > 300 keV, which is where the HB bounds stop. 

Also supernova explosions (SN) are used to put limits on axions. Stars with 
6-8 M mass or more reach the ultimate phase of the processing of nuclear fuel, 
creating an iron nucleus. Iron has the largest binding energy per nucleon and 
therefore cannot be efficiently burnt inside a star. It does not contribute to produce 
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the radiation pressure necessary to contrast the gravitational pull and to maintain 
the hydrostatic equilibrium. If the iron core reaches a critical mass, it collapses 
under its own weight. Once the nuclear density is reached the collapse stops and 
the bounce produces a shock wave that expels the outer layers in a core-collapse 
SN explosion. In the collapse, electrons are jammed inside protons forcing inverse 
/5-decays +p + — > n + u e : a neutron star forms and lots of neutrinos are created. 
The density of matter in a SN core is so high, that even neutrinos remain trapped 
and it takes some time before they can diffuse out [46]. After SN 1987A, 24 neutrino 
events were measured above the background in about 10 s. Their distribution in 
energy and time agrees well with the standard picture for type II SNe. A particle 
with a weaker matter interaction than the neutrino would provide a more efficient 
energy dissipation channel than the standard ones. If this is the case, the neutrino 
burst duration would have been shorter than what was measured. In such a high 
nucleon density environment, axions would be produced by virtue of their nucleon 
coupling in reactions like N + N—>N + N + a. In figure 1.2, the upper bound on 
the axion mass of SN 1987A, m a < 16 meV, thus /„>4x 10 8 GeV, comes from 
this burst duration argument. Lighter axions, and thus less coupled ones, would 
not be efficiently produced and thus they would have not significantly affected the 
timing of SN 1987A neutrino events [49]. On the other side, if axions couple much 
more strongly to matter, they could be trapped inside the core, and the neutrino 
burst would have suffered little or no modification. This is why the SN 1987A 
bound stops at low f a . However, in this case some axions are emitted and it would 
have been possible to detect them in the Kamiokande II experiment thanks to the 
nuclear reaction a + 16 -> 16 0* -> 16 + 7. The region labelled Kamioka 16 is 
excluded by the non-observation of these events [53]. The SN 1987A bounds are 
very uncertain and have to be taken with a grain of salt. They are particularly 
interesting if the axion has not a direct coupling with the electron, for the stronger 
limits coming from HB stars are not valid in this case. 

Most of the bounds just described are also valid in the ALP case, as they are 
directly constraining g a77 and thus g^. Moreover, those based on measurements 
of the coupling to electrons or nucleons can also be used to constrain ALPs, once 
they are translated to the effective ALP-fermion coupling mediated by photons. 
Because of their larger parameter space, ALPs have a richer phenomenology and 
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Figure 1.4: ALP limits. The blue regions are constrained by particle 
physics experiments. In green are the astrophysics related limits, in particular 
CAST+SUMICO are helioscope measurements. The laser experiments exclude 
the yellow region. Haloscopes have explored the brown region and in red are the 
bounds provided by cosmology. The dotted ellipses suggest the values of the pa- 
rameter space which should characterise a pseudoscalar particle in order to solve 
some astrophysical conundrums. The band of axion models includes values of C 7 
between 0.6 and 6. 



more constraining experiments and observations than the axion. Figure 1.4, where 
we collected the most relevant ALP bounds, provides a panoramic view of the ALP 
parameter space, in comparison with the black-dashed band which represents the 
axion models, whose central value is the KSVZ axion and the shading covers the 
interval 0.6 < C 7 < 6. 

The bounds related to cosmology are plotted in red and they will be explained 
in the next chapters, together with the brown constraints that come from the 
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search for ALP DM with ADMX and other haloscopes. In blue are the bounds 
coming from particle physics experiments, in green astrophysical observations and 
arguments. The yellow area is constrained by laser experiments dedicated to the 
search for ALPs. 

Particle physics experiments, as in the axion case, can set constraints only for 
quite strongly coupled particles and they pile up in the upper part of figure 1.4. 
They are very useful anyway, for they can test larger mass than astrophysics. 
In figure 1.4 we draw in blue the exclusion bounds due to T decay and positro- 
nium annihilation into invisible channels, labelled respectively T — > 7 + inv and 
e+e~ — > 7 + inv. The beam dump experiments performed in SLAC exclude the 
patch labelled Beam dump [54,55]. 

Astrophysics plays a central role in constraining the ALP case. The HB bound 
on c/ ar) , described before is directly applicable to ALP photon coupling, and it is 
plotted in figure 1.4 in light green and labelled HB. The SN 1987A gives some 
constraints too. The first one, labelled SN, is given by the duration of the mea- 
sured neutrino pulse, as in the axion case, and was derived in [30]. The second 
bound related to SN 1987A is labelled SN 7 burst and it is related to the trans- 
parency of the dense SN core to the ALP propagation. If an ALP exists in this 
region of the parameter space, it would be produced during the core collapse and 
it would subsequently escape from it. Then, the propagating ALPs can oscillate 
into high energy photons interacting with the galactic magnetic field, and finally 
detected on earth. Since no 7-ray pulse was measured in correspondence of SN 
1987A, the SN 7 burst region can be excluded [56]. However the SN bounds 
can be considered rather weak, as they rely on an insufficient understanding of 
the SN dynamics and, in the cases just discussed, of the ALP emission from a 
nuclear-density environment. 

The CAST helioscope constrains a large part of ALP parameter space. The 
Japanese experiment SUMICO gives also a bound in a small part of the parameter 
space not constrained by CAST [50]. Their limits are both plotted in dull green 
in figure 1.4 and labelled CAST+SUMICO. It is interesting to notice in this 
picture how it is actually difficult to reach the sensitivity to constrain the axion. 
In particular, just a corner of the bound penetrates the axion model stripe, which 
implies the slimness of the CAST bound in figure 1.2. 
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A light shining through the wall (LSW) experiment is the table-top version of 
an helioscope, in the sense that in this case the source is not the sun but a high 
intensity laser. The experiment consists of a pipe divided in two by an opaque 
barrier, the wall, in a production section and a detection one. In both of them a 
high vacuum has been created and an intense magnetic field imposed. The beam 
of photons emitted by the laser passes through the magnetic field before ending 
its run against the wall, converting some of the photons of the beam into ALPs 
via the Primakoff effect. Because of their feeble coupling, ALPs can cross the 
wall passing into the detection section, where they have to be converted back into 
photons interacting with the magnetic field before being finally detected [57]. The 
yellow patch labelled ALPS in figure 1.4 shows the constraints of the latest LSW, 
the Any Light Particle Search (ALPS) [58]. 

Measuring changes in polarization is a different approach to laser experiments 
[47,59]. Again, the laser beam passes through an intense magnetic field. The 
component of the vector potential parallel to the magnetic field effectively mixes 
with the pseudoscalar field and it is partially absorbed and retarded. This causes 
a light beam, linearly polarised at a given angle with the magnetic field, to rotate 
a small amount because of the absorptive process and to gain a small ellipticity 
due to the dispersion effect. Therefore, if pseudoscalar particles coupled to two 
photons exist, a magnetic field induces respectively dichroism and birefringence 
on the magnetised vacuum. Both effects are proportional to the square of the 
magnetic field and of the photon coupling g^. Moreover, the acquired ellipticity is 
proportional to m| [60] . The QED background processes are virtual pair creation 
for birefringence [61] and photon splitting for dichroism [62]. Both of them are very 
suppressed, especially the photon splitting, thus measuring significant magneto- 
optical properties of the vacuum would be a signal of the existence of particles 
coupled to two photons. The best constraints on these phenomena are provided 
by the Polarizzazione del Vuoto con LASer (PVLAS) experiment, whose results are 
shown in figure 1.4 in the yellow zone labelled PVLAS [63]. In particular, this 
refers to the birefringence measurements, for the dichroism ones are completely 
superseded by ALPS results. 

Finally, the MWD region is constrained by similar consideration, but powered 
by the high energy phenomena which sometimes characterise the astrophysical 
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environments. In the past decades many white dwarfs with very strong magnetic 
fields — up to 10 10 G — were discovered. We call these objects magnetised white 
dwarfs (MWD). The polarisation of light coming from MWD is strongly influenced 
by the magnetic field. A typical fraction of 5% of the light is circularly polarised, 
as required to electromagnetic waves propagating in a magnetised atmosphere, 
while the linearly polarised fraction is something less but of the same order, and 
could be explained by photon-ALP conversion. If ALPs with < 10~ 6 eV and 
94> ~ 10~ n GeV -1 exist, the linearly polarised fraction of light would be larger 
than 5% for MWDs with B = 10 9 G, and consequently this region has to be 
excluded [64]. The bound can improve by more than an order of magnitude for 
the coupling to photons — up to ~ 1CT 12 GeV -1 — if the data about the MWD 
with the strongest magnetic field, B ~ 10 10 G, are confirmed [64]. 

1.4 Where could axions and ALPs hide? 

Until now each axion and ALP search has been unfruitful, providing only exclusion 
bounds. However, besides the phenomenological need to have an axion in order 
to solve the strong CP problem, there are some observational hints about which 
regions of the parameter space could hide an axion or an ALP. We draw dashed 
ellipses in figures 1.2 and 1.4 to highlight the regions of parameter space suggested 
by these observations. 

Recently, the measurement of TeV photons from very far sources — namely 
some active galactic nuclei (AGN) — has puzzled astrophysicists [65,66]. Very high 
energy photons should inelastically scatter with the background light producing 
e + e _ pairs. They should be rapidly absorbed by the intergalactic medium and 
we should observe none of them if coming from very far sources. Of course, it 
could be that the spectrum of the source is much harder than we expect or that 
we overestimate the amount of background photons in the medium. However, it 
could also be that these TeV photons, interacting with the extragalactic magnetic 
field, oscillate into ALPs, then escape the absorption and finally, once converted 
back into the photon form, are detected by our telescopes [67-69]. To solve this 
transparency issue, this mechanism needs an ALP with very low mass, m$ ~ 
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10~ 12 -10~ 9 , and coupling greater than ~ 10~ 12 GeV -1 [56]. Of course, much 
about these considerations depends on the assumptions on the source and on the 
extragalactic magnetic field. Especially the strength of the latter has to be assumed 
very close to the upper limit to have an appreciable effect. To avoid this problem, 
one can assume the photon-ALP conversion to happen in the regions around the 
source and in our galaxy, where the magnetic fields are better known [70,71]. We 
have nevertheless to notice that if the constraint coming from MWD is solid, not 
too much space is left for ALPs to solve the transparency of the universe problem. 
But even if this is the case, the ALP-photon conversion affects the polarisation of 
distant astrophysical sources and, once some more information about the coherent 
components on the intergalactic magnetic field are obtained, it will be possible 
to extract some useful limits from AGNs [72, 73] and other celestial objects, like 
quasars [74]. 

Astrophysics provides another very interesting clue related to the evolution of 
white dwarf stars. It seems that if the axion has a direct coupling to electrons 
and a decay constant f a ~ 10 9 GeV, it provides an additional energy-loss channel 
that permits to obtain a cooling rate that better fits the white dwarf luminosity 
function than the standard one [75]. The selected mass range is in the meV range 
and g ail ~ 1CT 12 GeV -1 . The hadronic axion would also help in fitting the data, 
but in this case a stronger value for g ail is required to perturbatively produce an 
electron coupling of the required strength. 

Finally, axions and ALPs are also perfect dark matter candidates. We will 
devote a large section of chapter 2 to this topic. 

After the presentation made in this chapter, we will deal from now on with the 
cosmological bounds. Since many of the topics that will be treated involve both 
axions and ALPs, we will refer to them using the term pseudoscalars when they 
are on the same level. We also want to underline that in the following discussion 
g ail and g$ are phenomenologically equivalent. 



CHAPTER 2 



Establishing an axion or ALP relic 

population 



In order to obtain information about the axion and its ALP relatives from cos- 
mological considerations it is first of all necessary to understand if a primordial 
population of pseudoscalars can be established. Several mechanisms can achieve 
this task and they will be described in sections 2.1 and 2.2. Sections 2.3 and 2.4 deal 
with the ways a pseudoscalar population can disappear from the cosmic plasma. 
In the rest of the dissertation we will deal with a flat Robert son- Walker metric 



where R is the cosmic scale factor, which has length dimensions, t is the time 
coordinate, r is the dimensionless radial comoving coordinate, and (■&, if) are the 
dimensionless comoving angular coordinates. The cosmic scale factor is growing 
in time, representing the expansion of the universe. A useful way of measuring 
the expansion is through the redshift z, which measures the ratio between the 
wavelength A e of a light signal emitted at time t e and the wavelength of the 
same signal detected at time t d , and it is defined to be 



The redshift of a signal measured today is also a practical way to refer to cosmo- 
logical time scales. From the present temperature of the CMB, T = 2.73 K = 



ds 2 = dt 2 - R 2 (t) (dr 2 + rW + r 2 sin 2 tf 2 V) , 



(2.1) 





(2.2) 
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2.35 x 10~ 13 GeV, the temperature at a given redshift is easily obtained with the 
formula T = T (l + z) if we assume no heating of the thermal bath. 
The expansion rate is defined to be 

Hit) . || (2.3, 

where the dot stands for the time derivative. The Friedmann equation links H 
with the energy density of the universe p, 



H 2 = — -V , (2.4) 
3 m?' 



57T p 

the Planck mass being m P \ = 1.2211 x 10 19 GeV. In the radiation dominated 
universe, the energy density is 

P = ^g*(T)T\ (2.5) 
which depends on the temperature T, and the equation (2.4) becomes 

H-lMg^T) 1 ' 2 — . (2.6) 
m P i 

We will refer very often to this form for H. The T dependent quantity g* is the 
number of relativistic internal degrees of freedom, which is plotted in figure 2.1. 
Its definition is 

i=bosons ^ ' i=fermions ' 

where the indices % run over the bosons and fermions with temperature Tj and gi 
internal degrees of freedom, which are relativistic when the photon temperature 
is T. 

The critical energy density 

(£^)! = 105;!2 keV_ 
87T cur 

where H = 100 h kms _1 Mpc _1 is the present value of the expansion rate and 
h ~ 0.7 is its present-day normalized value, defines the energy density of a flat 
universe expanding at H rate. It can be used as a unit of measure for the energy 
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Figure 2.1: The number of relativistic internal degrees of freedom (solid) and 
the number of relativistic entropy degrees of freedom (dashed) as functions 
of temperature. 

density of the different constituents of the universe. Thus, defining the present 
ratios Q r = p r /p c , Vt m = p m / p c , and ^a = Pa/Pc respectively for radiation, matter 
and vacuum energy, the expansion rate can be conveniently expressed as a function 
of the redshift with 

h(z) = H oy /n r (i + z y + n m (i + zf + n A , (2.9) 

if there are no reactions converting one energy form into the other. 
The entropy density of the radiation dominated universe is 

a = %9*s{T)T* , (2.10) 
and g*s are the number of relativistic entropy degrees of freedom, 

9*3{T)= £ 9*(f)+\ E 4f) , (2.H) 

i=bosons ^ ' i=fermions 

which is also plotted in figure 2.1. For all the usual cosmological quantities, like g* 
and g*s, we follow the definitions of "The Early Universe" by Kolb and Turner [76]. 
We also define comoving quantities scaling out the effect of the universe expansion. 
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Therefore, comoving quantities change only in force of physical processes which 
are not the stretching of space-time. Some examples are the comoving entropy 
S = sR 3 , the comoving number density N = nR 3 and comoving momentum 



2.1 Production mechanisms I: non-thermal relics 



Soon after the invisible axion proposal, it was recognised that this new particle is 
a perfect DM candidate. The spontaneous breaking of the PQ-symmetry and the 
consequent phase-transition, which shifts the vacuum of the theory to the minimum 
of the potential (1.5), provides the axion with an efficient non-thermal production 
mechanism [77-79]. This realignment mechanism was later proposed also for the 
non-thermal productions of string axions [80-83], general ALPs [84] and hidden 
photons [85]. More recently, the realignment mechanism for hidden photons and 
ALPs has been reviewed in [3]. Here it is also shown that once produced, a 
population of very light non-thermal dark matter particles is extremely difficult to 
reabsorb in the primordial plasma. 

2.1.1 Realignment mechanism 

In the realignment mechanism a field, which in the early universe can take a 
random initial state, rolls down towards the minimum of the potential. Once it 
has reached the bottom, it overshoots the minimum and starts to oscillate around 
it. If the quanta of the field are cosmologically stable, these oscillations behave as 
a cold dark matter fluid. Their energy density in fact is diluted by the expansion 
of the universe as p oc R~ 3 . 

The simplest example is that of a scalar field of mass with Lagrangian 



where Cj encodes interactions of the scalar field with itself and the rest of particles 
in the primordial bath. In general, the mass receives thermal corrections from Cj 
which might be crucial, thus = m^it) should be understood. In each causally 
connected patch of the universe the scalar field has an initial value, 0j. If inflation 



K = kR. 



C = - \ 



(2.12) 



2.1 Production mechanisms I: non-thermal relics 



29 



1.0 

0.5 
0.0 
-0.5 

0.5 1 5 10 

Log 10 (t/fi) 

Figure 2.2: Numerical solution of equation (2.13) (solid) and the evolution of 
log 10 (m^/H ) (dashed) for the values used in the equation. 

already happened, is uniformly equal to 0j in the whole observable universe. The 
equation of motion in the expanding universe for the homogeneous component of 
— called the zero mode — is obtained neglecting the gradient effects, 

+ 3#0 + mj0 = . (2.13) 

Its solution can be separated into two regimes. In a first epoch, 3H ^> m^, so is 
an overdamped oscillator and gets frozen, = 0. At a later time, ti, characterized 
by 3H(ti) = m^ti) = m 1; the damping becomes undercritical and the field can 
roll down the potential and starts to oscillate. During this epoch, the mass term 
is the leading scale in the equation and the solution can be found in the WKB 
approximation, 




(2.14) 



where 0(ti) ~ 0j, since up to t\ the field evolution is frozen. Note that to obtain 
this solution only the definition of H = R/R has been taken into account, not 
its actual time dependence, and so it is valid for radiation, matter, and vacuum 
energy dominated phases of the universe and their transitions. Figure 2.2 shows a 
numerical solution of equation (2.13). 
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The approximate solution (2.14) corresponds to fast oscillations with a slow 
amplitude decay. Defining this amplitude A(t) = ^(mi-Rf/m^i? 3 ) 1 / 2 and the 
phase a(t) = f* m^tydt, the energy density of the scalar field is 

P4> = + l -ml<p 2 = l -mlA 2 + ... , (2.15) 

where the dots stand for terms involving derivatives of A, which by assumption 
are much smaller than m^, because ^> H in this regime. The pressure is then 
1-1 1 

P<f> = ~ 9 m |^ 2 = _ 9 m |^ 2 cos ( 2a ) _ AAm^ sin (2a) + A 2 cos 2 (a) . (2.16) 

When the field just starts to oscillate the equation of state is a non-trivial and 
strongly time dependent function. However, at much later times, t t±, the 
oscillations in the pressure occur at time scales much faster than the cos- 

mological evolution. We can therefore take an average over these oscillations, and 
the pressure is then 

(^) = a 2 cos 2 (a)) = ^ 2 . (2.17) 
At leading order in A/ (Am), the equation of state is just 

w = (p)/(p) ~ , (2.18) 

which is exactly that of non-relativistic matter. 

It follows from equation (2.14) that the energy density in a comoving volume, 
p<f>R 3 , is not conserved if the scalar mass changes in time. However, the quantity 

P<t> R3 = 1^ D3^2 

2 

is constant, and can be interpreted as a comoving number of non-relativistic quanta 
of mass 1710. We can use this conservation to compute the energy density today to 
obtain 

PSo) = ™o^| * , (2.20) 

where quantities with a 0-subscript are evaluated at present time. 

Through the conservation of comoving entropy S, the dilution factor (Ri/Kf 
in equation (2.20) can be expressed as 

RA 3 g* s (T)T 3 



N <t> = ^ = >!i??0 2 , (2.19) 



RJ "^(W (2 ' 21) 
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Figure 2.3: The function T(T), used in equations (2.22). 

Using the expression for the Hubble constant in the radiation dominated era, 
equation (2.6), and the definition of T\, 3H(Ti) = mi, T\ can be written in terms 
of mi and the Planck mass. In this way, equation (2.20) for the matter density from 
the misalignement mechanism today can then be expressed in a more quantitative 

way as 

— (°-S)^(lo4v) 2 ^). W 

(ira) 2 -^ - (2 - 22b) 

where T(Ti) = (g i) ,(Ti)/3.36)i(g^s(Ti)/3.91)~ 1 is a smooth function ranging from 
1 to ~ 0.3 in the interval T\ e (T ,200GeV), and it is plotted in figure 2.3. The 
quantity VL^ = p^o/pc is the fraction of energy density of the universe in ALPs. 
The abundance is most sensitive to the initial amplitude of the oscillations, being 
proportional to (f)f, and to a lesser degree to the present mass value m . The factor 
1 / ^fm{ reflects the damping of the oscillations in the expanding universe: the later 
the oscillations start, i.e. the smaller Ti and therefore Hi and mi, the less damped 
they are for a given m . 

If we compare the above estimate with the DM density measured by WMAP 




^0 ~ [0.016 h ) W-rrW — 
' V e V y mi 
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and other large scale structure probes [86], 

keV 



Pcdm = 1-17 — - , Qcdm = 0.11 hT 1 , (2.23) 



cm 



it seems that only very large values of fa can provide a p^ that accounts for all 
the dark matter. However, a relatively small fa could be compensated by a small 
mi < m . 

It is now necessary to provide a model for the thermal mass of the pseudoscalar 
particle to compare the estimate (2.22a) with the observed abundance (2.23). For 
the potential (1.25), the ALP satisfies the equation of motion (2.13) as long as the 
angle if = (/)//</> is small. The inaccuracy of the quadratic approximation can be 
anyway cured by an additional correction factor to (2.22a). This is normally an 
(9(1) factor except if we fine tune the initial condition to ip — ir. 

Figure 2.4 represents the region of the parameter space of ALPs in which the 
realignament mechanism could provide a sufficient production of DM. The allowed 
regions of ALP dark matter in the m^-g^ plane in figure 2.4 can be obtained using 

4 = Vi ^L (2.24) 

with ifii = fa/ fa, the initial misalignment angle, whose range is restricted to values 
between — n and +it. The model dependent factor C 7 will from now on be taken 
to be unity for presentation purposes, but the reader should keep in mind that 
in principle it can assume a very different value according to the model it comes 
from. To tune the right DM abundance we have used the a priori unknown value 
of ifi. The upper bound on reflects the fact that ip^ cannot be larger than 
7r, and thus assumes ipi ~ ir. Moving to lower values of requires inflation 
happening after SSB in order to have a homogeneous small value of ipi which is 
increasingly fine-tuned to zero to avoid over-abundant DM production. In this 
sense the values closest to the boundary, corresponding to the largest values of 
the photon coupling, can be considered the most natural ones. We will further 
comment on the fine-tuning issue later on. 

In the simplest realization of an ALP model, the mass receives no thermal 
corrections, thus is constant throughout the universe expansion and the DM 
yield can be inferred directly from equation (2.22a) using m\ = m . In figure 2.4 
it is the dark red region labelled m = mi. 
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Figure 2.4: Axion and ALP DM from the realignment mechanism. The axion 
model band is plotted as in figure 1.4. The regions of the ALP parameter space 
in which the realignment mechanism provides sufficient — or too large — DM 
abundance are plotted in different red gradations for different thermal mass mod- 
els. The pink region fulfils the minimum requirement for non-thermal produced 
ALP to behave like CDM. The grey region cover the part of parameter space in 
which the ALP lifetime, which is provided by equation (2.36), is shorter than the 
age of the universe. The turquoise regions have been tested by haloscopes and 
no ALPs have been found. As a reference, the HB upper bound is plotted as a 
dotted line. 
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If the global symmetry associated to the ALP is anomalous, then the mass is 
provided by the instantonic potential just as the rf or the axion acquire their mass 
thanks to QCD instantons. The ALP needs in principle another unbroken SU(N) 
group, which condenses at a scale A. If T is the temperature of the SU(N) sector, 
then the thermal ALP mass is 

= m for T < A , 
m„~{ (2.25) 

m (ff forT>A. 

It makes sense to naively assume A' ~ A" ~ A, while the precise relation among 
these scales depends on the details of the model. At temperatures larger than 
A, electric-screening damps long range correlations in the plasma and thus the 
instantonic configurations, resulting in a decrease of the ALP mass. In specific 
models the exponent /3 can be obtained for instance from instanton calculations, 
but here it is left as a free parameter. Assuming the onset of ALP coherent 
oscillations to happen in the mass suppression regime, it is easy to obtain an 
expression for mo/mi which is the expected enhancement in the DM abundance, 

(^)*(3xl.66 v ^T)* (2.26, 

and the factor that controls the enhancement is 

y/m m P1 _ A' /m^T 




A" A"V n 



These models can provide only a moderate enhancement of the DM density with 
respect to the constant m<p case. The gained regions for the ALP DM case for values 
of P — 1, 3, 5, 7, 9 can be seen in figure 2.4 from bottom to top — the lowermost 
region m\ = m corresponds, of course, to j3 — 0. Actually, even considering 
unrealistically huge values of (5 does not help much, as can be seen from the 
asymptotic approach of the highest (5 cases. This is reflected by the finite limit of 
equation (2.26) when f3 — >■ oo, but it follows from its definition, equation (2.25). 
In the (3 — > oo limit, is a step function of temperature, oc 6(T — A), and 
the relation m = A 2 //^ determines A from m and fy. Thus each point in the 
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m <P~9<t> parameter space has an implicit maximum DM abundance, independent 
of p. The crucial assumption that leads to these conclusions is that A' ~ A, 
because it does not allow to consider arbitrary small values for A for a given mass. 
Therefore, models in which A' ^> A imply generically higher DM abundance and 
require smaller initial amplitudes fa. Unfortunately, at the moment we can not 
provide a fully motivated example. 

Finally, the minimal requirement for an ALP to behave as CDM is that at latest 
at matter-radiation equality, at a temperature T eq ~ 1.3 eV, the mass attains its 
current value m and therefore the DM density starts to scale truly as 1/R 3 . In 
particular, at this point the field should already have started to oscillate. This 
corresponds to a lower limit on mi, mi > 3H(T cq ) = 1.8 x 1CT 27 eV, which implies 
an upper bound on p^ i0 , 

In other words, if we want these particles to be the DM, we need 

(m /eV)(&/54TeV) 2 > 1 , (2.29) 

giving a constraint on the required initial field value as a function of the mass today, 
which is plotted in figure 2.4 in pink. Again, we are not able to provide a well 
motivated model for an ALP that could perform this late realignment production 
and this limit has to be considered just as the broadest region of ALP parameter 
space that in principle could provide the right amount of CDM. 

The dark matter generated by the realignment mechanism has interesting prop- 
erties beyond those of cold dark matter. At the time of their production, particles 
from the realignment mechanism are semi-relativistic. Their momenta are of the 
order of the Hubble constant p ~ H 1 < Accordingly we have today — out- 
side of gravitational wells — a velocity distribution with a very narrow width of 
roughly, 




(2.30) 



Combined with the high number density of particles, n^ = N^/Rq = Pcdm/^o, 
this narrow distribution typically leads to very high occupation numbers for each 
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quantum state, 

Occupation ~ ^ ^ ~ 10 ^ — j (^-J , (2.31) 

where we used R /Ri ~ Ti/T ~ a/^i^pi/T . If the ALP self-interactions are 
strong enough to achieve thermalisation and ALP-number-conserving, as argued 
in references [87, 88] for the case of axions, this high occupation number leads 
to the formation of a Bose-Einstein condensate. This could imprint interesting 
signatures in cosmological observations, like galactic caustic rings [87-90]. 

The axion realignment production is a particular case of what is described 
above. The A parameter in this case is Aqcd, and the phase transition that 
turns on the axion mass is the QCD phase transition. The axion thermal mass 
was first calculated from the dilute instanton gas approximation [91] and more 
recently using the instanton liquid one [92]. The theoretical uncertainties in these 
calculations make the parameters for the thermal axion mass slightly different. For 
example, a simple and recent approximation in the dilute gas approximation, that 
also agrees very well at high temperature with the liquid instanton one, is [92] 

A 2 /A\ 3 ' 34 

m a (T) = 4.1 xl(T 4 — (^-J , (2.32) 

for A = 400 MeV. This approximation is valid above T ~ 0.1 GeV, since below 
this value the axion mass approaches the zero-temperature value m a given by 
equation (1.6). The exact shape of the thermal axion mass determines ti, thus it 
affects the m Q /mi ratio in equation (2.22a). 

2.1.2 Before or after inflation? 

We do not know when cosmic inflation occurred and if it took place before or after 
the spontaneous breaking of the PQ-symmetry. This is one of the crucial points 
in calculating the axion energy density from the realignment mechanism. If the 
initial (pi was set after inflation, we need an average of y? 2 over all the universe to 
calculate the amount of axions produced through equation (2.22a), for the current 
universe includes many patches that were initially causally disconnected. However, 
this is not possible if the universe inflated at T < f a : in this case the resulting 
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axion energy density p a remains </? 2 dependent. Assuming a flat prior on the value 
of the initial misalignment angle, ifi ~ 0(1) is the most likely outcome. If this 
is the case, it is possible to state that the axion mass should be m a > 6 /xeV, 
and thus f a < 10 12 GeV, to avoid DM over-production. This is the traditionally 
quoted lower bound, first obtained in 1982 [77-79]. It is represented in figure 1.2 
with the Cold DM label. Successive estimates, with a more refined treatment of 
the average of the squared initial angle or of the thermal axion mass, give more 
severe bounds but of the same order. For the realignment contribution, the recent 
reference [92] gives a limit which is stronger by one order of magnitude. The 
traditional f a < 10 12 GeV bound is therefore the most conservative one. 

It seems that this cold DM bound could be easily avoided in the case inflation 
happened after the spontaneous symmetry breaking, if we accept a fine-tuned very 
small value for </?j. In figure 2.5 are plotted the isocontours of ifi that provide the 
whole cold DM density, if in equation (2.22b) we take m x = m and T{T) = 1. 
However, ifi has a minimum value due to the quantum fluctuations generated in 
the pseudoscalar field during inflation. These unavoidable inhomogeneities in <fi 
are of order Sifii ~ iff/ (2-71"/^), where Hj is the value of the expansion rate at 
the end of inflation. This precludes fine-tuning of the universe average of </>j below 
if//27r, and sets a minimum DM abundance for a specified value of Hj. Since the <fr 
field is effectively massless during inflation in this scenario, these inhomogeneities 
correspond to isocurvature perturbations of the gravitational potential. This has 
been discussed extensively in the literature in the context of axions and of string 
axions, see for example [81,93,94]. 

The WMAP7 observations of the primordial density fluctuations set very strin- 
gent constraints on isocurvature perturbations from which one can obtain an upper 
bound on Hj, assuming a given f^. WMAP measures [86] 



where (IS* 2 !) ~ is the isocurvature power spectrum, and (|-R 2 |) the adiabatic 
one, which is generated by the inflaton or by other fields. At the pivot scale 
k = 0.002 Mpc" 1 , WMAP finds (|i? 2 |) = 2.42 x 10" 9 , which gives the bound 



a = 




< 0.077 , 



(2.33) 



Hj < 4 x 10-Vi/* . 



(2.34) 
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Figure 2.5: The anthropic window. The solid lines represent the isocontours 
of the values of (fi that provide the whole cold DM density, according to equa- 
tion (2.22b), taking mi = tuq and F{T) = 1. Dashed lines give the isocontours of 
the maximum value of Hi allowed by the limit (2.34), calculated using the same 
value of ifi as before. As a reference, the KSVZ axion parameters are plotted as 
a dotted line. In the grey region ALPs are cosmologically unstable. 



The lower bound for the reheating temperature is T RH > 4 MeV [95]. Linking 
naively the two quantities gives H I « T^ H /m P1 > 1.3 x 1CT 24 GeV. In the light 
of equation (2.34), it is rather a weak bound by itself. However, in the most 
common inflation models, Hi is preferred to be much larger, to allow mechanisms 
that require very high initial temperatures of the universe — like leptogenesis, 
which needs T > 10 9 GeV [96]. The isocontours of the maximum Hi allowed by 
limit (2.34) are plotted in figure 2.5, where the value of ^ that provides the whole 
DM density is used. A given value of Hi excludes the region of the parameter space 
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that lies on the right of its isocontour. Then, along the Hj isocontour selects 
the tfi that is provided by primordial fluctuations, while on the left half-plane the 
initial misalignment allowed is larger and thus less fine-tuned. 

In the case of the axion, the values of f a that satisfy the constraint (2.34) for 
large H I are above the GUT scale, which seems very reasonable from the model 
building point of view: they form the so-called anthrophic window 1 . Considering 
that the axion was introduced to avoid fine-tuning to solve the strong-CP problem, 
the legitimacy of the choice of a very small ipi to justify an axion scale linked to 
Agut or the string scale is questionable [97]. However, the anthropic point of view 
can find a vindication regarding the DM density and other significant physical 
quantities simply as observational data, from which it is possible to extract the 
value of the initial misalignment with an a posteriori analysis [98]. This procedure 
is legitimate because of the intrinsic randomness of ipi, and converts the fine-tuning 
issue into an observational problem. 

On the other hand, if the SSB took place after inflation the fa picked up random 
values in different causally disconnected regions of the universe and the anthropic 
window loses its significance. If this is the case, also the gradient effects should be 
taken into account. At larger scales, the DM density averages to a constant value 
corresponding to (faf) ~ 7r 2 /|/3, bearing the mentioned (9(1) correction due to 
the non-harmonic behaviour of large initial phases. The initial size of the domains 
of different fa cannot be larger than 



Non-linear effects, due to the attractive self-interaction caused by higher order 
terms in the expansion of the potential (1.25), drive the overabundances to form 
peculiar DM clumps that are called miniclusters [99-102]. These act like seeds 
enhancing the successive gravitational clumping that leads to structure formation. 
The minicluster mass is set by the dark matter mass inside the Hubble horizon 
d H = H^ 1 when the self-interaction freezes-out, i.e. M mc ~ p ( p(T\)d H (T\) 3 for 

1 The adjective anthropic refers to the fact that the existence of the human species requires 
the universe to fulfil some conditions about its age and composition, which too abundant axion 
DM would not meet. 




(2.35) 



H SSB f}y/g*(U) 
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the freeze-out temperature T\. Long-range interactions will be exponentially sup- 
pressed at distances longer than l/m^ so we can expect T\ to be of the order T l5 
with at most a logarithmic dependence on other parameters. This is indeed the case 
for QCD axions, for which the miniclustering is quenched soon after the QCD phase 
transition that turns on the potential (1.25) [87] giving M mc ~ 1O~ 12 M , where 
M Q = 1.989 x 10 33 g is the solar mass, and a radius R mc ~ 10 11 cm [103]. In the 
case of ALPs, M mc can be larger if the mass is lighter. The authors of [102] pointed 
out that the present data on the CDM power spectrum constrain M mc < 4 x 1O 3 M 
which translates into a lower bound in temperature T\ > 2 x 10~ 5 GeV and in the 
ALP mass > H(T = 2 x 10~ 5 GeV) ~ 10~ 20 eV. If some of these miniclusters 
survive the tidal disruption during structure formation they should be observable 
in forthcoming lensing experiments [102,103]. 

2.1.3 Contribution of topological defects 

During the spontaneous symmetry breaking of a global symmetry topological de- 
fects such as cosmic strings and domain walls are formed [104]. Strings form 
because of the breaking of a U(l) symmetry. They have a thickness 5 ~ 1//^ and 
typical sizes of the order of the horizon, L ~ t. As strings enter into the horizon 
they can rapidly reconnect, form loops and decay into PNGBs and we have to 
consider also their contribution to the axion and ALP energy density. 

Axions resulting from string decay are known to contribute significantly to 
their cold DM density, but the exact amount is subject to a long-standing con- 
troversy [76]. The debate is focused around the axion emission spectrum. Some 
authors argue that the decay proceeds very fast, with an emission spectrum 1/k 
with high and low energy cutoff of order respectively 1/5 and l/L. In this case the 
contribution to is of the same order as the realignment mechanism [105, 106]. 
Others put forward that the string decays happen after many oscillations, with a 
radiation spectrum peaked around 2n/L, which enhances the contribution to cold 
DM by a multiplicative factor of log(L/5) ~ log(/ /m ) ~ 0(60) [107-113]. 

Once the axion potential builds up at the QCD phase transition, also domain 
walls build up. If the axion potential has only one minimum the domain walls 
can still efficiently decay into axions — since on their two sides there is the same 
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vacuum — leading to a third axion population. If different exactly degenerate 
vacua exist 2 the domain walls are persistent. The energy stored in a domain 
wall is huge, and can very easily run in conflict with observations. An option to 
overcome this problem requires to assume a small explicit breaking of the Peccei- 
Quinn symmetry, which breaks the degeneracy and select one among the possible 
vacua, otherwise it has to be Af = 1 [114]. It has been recently pointed out that 
the axions emitted by decaying domain walls constitute the dominant population 
in the slow string decay scenario [115]. Summing all the contributions to axion 
non-thermal DM in this scenario would constrain the axion decay constant to be 
f a < (1.2-2.3) x 10 10 GeV, two orders of magnitude lower than the traditional 
axion CDM bound [115]. We plotted this bound in figure 1.2, only enclosed by a 
line because of its uncertainty: in the case the SSB would have happened before 
inflation, the topological defects, being diluted away, would provide a negligible 
contribution and thus the traditional value is again the one to be considered. The 
same type of behaviour is expected for ALPs with characteristics similar to those 
of the axion, i.e. ALPs whose mass is generated at a late phase transition due to a 
hidden sector which becomes strongly interacting. In this case we should keep in 
mind the controversy of the string decay spectrum and assume an uncertainty of 
order log( f^/m^) in the DM abundance. The domain wall problem can in principle 
be solved by strong enough explicit breaking if Af > 1. Again, the contribution to 
ALP energy density from wall decay should mimic the axion case in the Af = 1 
case. 

2.1.4 Detection of pseudoscalar cold dark matter 

The axion and the ALP are optimal cold DM matter candidates, but only if they 
are cosmologically stable, which depends on their mass and couplings. The two- 
photon coupling provides a decay channel, and the consequent lifetime is 



Of course, if the pseudoscalar is massive enough other decay channels are kinema- 
tically allowed. In figure 2.4 the grey area hides the part of the parameter space in 
2 From the definition (1.14) of f a , the axion potential (1.5) has exactly N equivalent vacua. 




(2.36) 
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which the ALP and the axion have a lifetime shorter than the age of the universe, 
which is approximately 10 17 s. In the rest of the dissertation we will very fre- 
quently deal with r, probably the most important parameter for our cosmological 
considerations. 

Axion and ALP cold DM is currently searched with resonant cavities called 
haloscopes [29]. These are high quality resonant microwave cavities permeated by 
a strong and static magnetic field. The cavity is tuned to convert the pseudoscalars 
in the DM halo into photons exploiting the two-photon coupling. The signal is a 
monochromatic peak at the frequency of the mass of the particle, broadened by 
the virial distribution of particle velocity in the galactic potential [22,116]. 

A haloscope provides a power of pseudoscalars converted to photons equal to 

W = glVB 2 ^Cmm(Q c ,Q*) , (2.37) 

where B is the magnetic field, V the cavity volume, C the form factor for the cavity 
mode, and Q c and are respectively the quality factors for the cavity and for 
the pseudoscalar signal, which is basically the ratio of the energy of the signal over 
its spread in energy [116]. A technological challenge in building such apparatus is 
to balance the quality factor of the cavity. For the moment, the highest Q c that 
can be obtained is still about an order of magnitude less than quality factor of the 
expected signal. Of course, a very high magnetic field is preferable. The geometry 
of the cavity selects the range of masses to which the instrument is sensitive and 
also influences the power (2.37) through the volume factor. The design of the 
device has therefore to find a compromise between the two, considering that the 
geometry of the cavity has to be optimised to measure the most probable range for 
pseudoscalar cold DM, which in the axion case means very suppressed couplings. 
The cavity moreover has to be tunable, in order to span several values of m^. To 
explore a broad mass range takes a long time, for just one mass value per time 
can be tested. Finally, since the signal provided by pseudoscalar DM would be 
rather feeble, many techniques to suppress the noise in the apparatus are applied, 
in order to enhance the signal to noise ratio. 

The first haloscope was the Rochester-Brookheaven-Fermilab detector [117], 
followed by the University of Florida experiment [118], but both of them had not 
enough sensitivity to test the axion. Only a few years later ADMX was almost able 
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to reach the sensitivity to test the axion as DM in the 1.9-3.5 /xeV range [42]. ALP 
cold DM was tested and not found in the turquoise regions labelled Haloscope 
excluded in figure 2.4. In this mass range, for > 1CT 10 GeV -1 this is not 
surprising because we do not expect any ALP condensate. The haloscope tested 
regions are also plotted in figure 1.4 and labelled Haloscopes. A new experiment, 
called ADMX-HF, is planned to operate at Yale University: the plan is to scan the 
ALP parameter space in the range around 10-100 fieV, with enough sensitivity to 
reach the axion band [45]. 

2.2 Production mechanisms II: thermal relics 

We learned in the previous section that the higher the ALP/axion decay constant, 
the more abundant would be the final yield, it being proportional to (fif oc /|, see 
equation (2.22a). If however is not high enough to guarantee a sufficient non- 
thermal production, the ALP or the axion can be nevertheless a DM component 
as it would be more strongly coupled and more prone to thermal production. 

A population of particles with g internal degrees of freedom — which counts 
the polarisation and the particle-antiparticle states — and mass m in thermal 
equilibrium has number and energy densities and pressure given respectively by 



if the temperature is T. According to the fermionic or bosonic nature of the 
particle, the sign + or — holds in the previous formulae. In the case of negligible 
chemical potential, the previous equations in the relativistic limit T ^> m for 
bosons are 




(2.38) 



(2.39) 



(2.40) 




(2.41) 



(2.42) 
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(2.43) 



while for fermions are 



n, 



Peq 



■cq 




(2.44) 



(2.45) 
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(2.46) 



Peq 



In the non-relativistic limit, T<m, the quantum statistics makes no difference 
and we have 



The dynamics of a particle distribution is governed by the Boltzmann equation. 
The Boltzmann equation links the total time derivative of a distribution on the left 
hand side to the microphysics of particle interactions that lies on the right hand 
side. The effects of the universe expansion are included on the left hand side. If <7j 
is the cross-section for the scattering of the particle under examination with the 
i-species in the bath, the Boltzmann equation for number density can be written 
as 



assuming that the interaction products rapidly thermalise, which allows us to write 
the term n ca n l e(l on the right-hand side [119]. The thermal average is obtained 
integrating in momentum the matrix element for the interaction process multiplied 
by the distribution functions of the particles involved. 

Our aim is to calculate the rate at which an axion or ALP population can 
arise from the scattering of SM particles in the thermal bath during the early 
phases of the universe. The axion and ALP self-interaction can be neglected, 
being suppressed by f^ 4 . Then, we can consider only the interactions with particle 
species in thermal equilibrium. We will specify better which interactions with 




(2.47) 



p cq = Tn . 



p eq = mn , 



(2.48) 
(2.49) 




(2.50) 
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which particles in the following. Under this assumption, ri 1 is always n* q and the 
Boltzmann equation simplifies as 



n 



3Hn = -Y (n - n cq ) , (2.51) 



where T = £V Tj is the total interaction rate and the i-th partial one is defined to 
be Ti = (tjiv)n l eq . From equation (2.51) we infer that if T is much bigger than the 
expansion rate H the number density is forced to follow the equilibrium one. If 
the right hand side becomes negligible, n is only diluted by the expansion of the 
universe, and changes through the factor provided by equation (2.21) 

n (*)= n (**)(^) 3 • ^ 

We can consider Y/H > 1 as a rule of thumb to define the thermal equilibrium 
condition, while if Y/H < 1 the species under examination is decoupled. Both Y 
and H are temperature dependent. When a particle species loses thermal contact 
with the bath because all its interactions freeze out, it is said to decouple, and 
this moment is characterised by the temperature Tf Q which is the solution of the 
equation H(Tf ) = r(Tf Q ). At any later time between the freezing out and the 
decay, the pseudoscalar density can be obtained in first approximation thanks to 
equation (2.52) and the dilution factor (2.21) 

MT) = n c ;(T fo ) f . (2.53) 

Assuming the particle content of the standard model, there is a minimum pseu- 
doscalar yield which is given by the value of g*{T io > E EW ) 

Our simple criterion about Y/H is sufficient to give a first approximation of the 
axion thermal history according to the axion scale /„. The key events take place 
during the radiation domination era, and therefore H can be expressed through 
equation (2.6). It is conceptually convenient to divide the thermal history of the 
universe into three epochs, classifying them through the main axion interaction at 
work. 
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fa>T > AqcD 

We learned in the previous section that when the universe temperature falls 
below f a the PQ-symmetry is spontaneously broken and the axion, still mass- 
less, pops up. If f a is not too high, the axion coloured interactions depicted in 
figure 2.6 are efficient. The interaction rate is T ~ {a^/ fa)n co \ cx (a^//f)T 3 , 
where n co \ is the number density of coloured particles. A more precise esti- 
mate of T is given in [120]. The ratio T/H cx T. Thus at high T the axions 
are in equilibrium — if the required temperature was ever achieved — and 
they decouple when T becomes too small. To decouple during this phase, 
the axion must have a decay constant in the range f a ~ 10 8 -10 10 GeV, thus 
a mass around the 1-10 meV. If this is the case, the fraction of the universe 
energy density in thermal axions is negligible, Vt a h 2 ~ 10 _6 -10 -4 [120]. 



o 
o 




Figure 2.6: Relevant axion interactions for T > Aqcd- 



• Aqcd >T>m n 

At T ~ Aqcd gluons and quarks confine. At this stage, if not yet decoupled, 
axions can still interact with pions thanks to the effective Lagrangian (1.23) 
(figure 2.7), with T ~ T 5 /(f v f a ) 2 [121]. If the axion decay constant is /„ ~ 
10 4 -10 7 GeV, which in terms of mass is m a ~ 0.01-100 eV, axions decouple 
at this stage. The axion yield in this case is larger than in the previous one, 
and in principle it can be abundant enough to influence the evolution of the 
universe. The formation of LSSs like those we observe today requires the 
percentage of relativistic matter to be not too high, otherwise free streaming 
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light particles would prevent the gravitational binding of the smallest clumps. 
Together with the CMB measurements, LSS constrains the axion mass m a 
to be smaller 0.7 eV and thus f a > 9 x 10 6 GeV [121, 122]. In figure 1.2 
this upper bound is labelled as Hot DM. This bounds applies to a small 
range of axion masses, for its validity ceases if the axion is cosmologically 
unstable. In this case the cosmological effects of the axion needs other tools 
to be analysed, which will be discussed in the following chapters. 




Figure 2.7: Axion-pion interaction relevant for Aqcd > T > m- 



Once the pions become non-relativistic, their abundance is exponentially 
suppressed, as it can be seen from equation (2.47), and the axion-pion in- 
teraction freezes out. Axions are only left to interact electromagnetically 
with photons and charged particles. The interaction of axions during this 
phase is less efficient in keeping them in equilibrium, not only because the 
electromagnetic force is weaker than the strong one, but also because the in- 
teracting species are just a few and very diluted by the expansion. Figure 2.8 
shows how a charged particle interacts with an axion via the Primakoff ef- 
fect, which remains active until the e + e~ annihilation at T ~ m e . For this 
interaction, T a77 ~ ag^ 7J n e . The axion decay is also an electromagnetic 
interaction, with r = 1/r. In particular, this process is characterised by a 
freezing in temperature, since r has no T dependence and thus the ratio Y/H 
grows with time. Axions with a mass larger than (9(10) keV never freeze out, 
because for them H ~ 1/r before the Primakoff interaction with electrons 
decouples [1]. In this case, the thermal equilibrium condition implies their 
disappearance from the bath when they becomes non-relativistic. 
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Figure 2.8: Relevant axion interaction for T < m 



Of course, if the axion has direct coupling to other SM particles, like the 
electron as in figure 2.9, the final yield would be higher, because, according to 
equation (2.51), a higher total interaction rate V means a later decoupling. In the 
ALP case, we are mainly interested in the interactions involving the two-photon 
coupling (1.24). Additional interactions would delay the decoupling and increase 
the final abundance. 



Figure 2.9: Additional axion-electron interaction, which is relevant in the DFSZ 
model for T < m^. 

As already mentioned, the cosmological stability is the second key factor to 
obtain limits on the parameters of pseudoscalars from cosmological observables. 
If the axion mass is larger than 20 eV, we can see from equation (2.36) that its 
lifetime is shorter than the age of the universe. Considering our schematic repre- 
sentation of axion thermal production as a function of f a , it seems worthwhile to 
give some more details about the last processes, for they determine the abundance 
of cosmologically unstable axions. Moreover, the generalisation to the ALP case 
is straightforward. 

The rate of Primakoff production due to scattering on relativistic electrons was 




2.2 Production mechanisms II: thermal relics 



49 



computed in [123] to be 



2 

P 12 



log I ^— ) + 0.8194 

mi 



(2.55) 



where m 7 2 = 2naT 2 /3 is the squared plasmon mass in a relativistic electron- 
positron plasma and T the temperature. Taking the rate to be proportional to the 
number density of electrons, n e = 3((3)T 3 /tt 2 , we can generalise to a multicom- 
ponent plasma, using instead the effective number density of charged particles 

n q = lEM 2 n t = ^-g q (T)T 3 . (2.56) 

i 

The parameter g q {T) represents the effective number of relativistic charged degrees 
of freedom. The plasmon mass has also to be corrected by a factor m 7 oc g l J 2 . 
The Primakoff interaction rate becomes 



2 2 

-p OL ga-yy 7T Tl q 

1 p ~ 



T 2 , 

^ (2-57) 



log ( — J + 0.8194 



12 3C(3) 

The freezing out temperature for the Primakoff interaction is easily calculated 
from Fp/H = 1 , which gives 

r t0 ^ " 2 ^^ 12 3^( 10 " 9GeV "') 2 GeV. (2.58) 
ag ajJ 2 m Pl g q g q \ g a ~ n J 

In the ALP case, for values of < 2 x 10~ 9 GeV -1 , interactions freeze out at 
temperatures above the electroweak scale, where the particle content of the plasma 
is somewhat speculative. For instance in the minimal supersymmetric standard 
model scenario, above the supersymmetry breaking energy scale we have g* = 
228.75, while the SM alone provides only 106.75 relativistic degrees of freedom. 
For < 10~ 17 GcV -1 we require a freeze out temperature above the Planck scale, 
which is most likely meaningless. 

Cosmologically stable ALPs must not exceed the measured abundance of DM, 
therefore 

n h 2 = ^h 2 = H^±h 2 < n DM h 2 = o.n . (2.59) 

Pc Pc 

Using (2.53), we see that ALPs with a mass = 154 eV would account for all 
the dark matter of the universe, thus larger masses are excluded. The discovery of 
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new degrees of freedom (dof) above the EW scale would relax this bound, which 
is linearly sensitive to g*(Tf ), by a factor (106.75+ new dof)/106.75. This means 
that (9(100) of them are needed for a sizeable change. 

2.3 Absorption of overabundant pseudoscalars 

The prototypical Boltzmann equation (2.51) shows that if an interaction rate is 
effective, the particle number density tends to its equilibrium value. Therefore, in 
case that the pseudoscalars are more abundant than those required by the equi- 
librium condition because of the realignment production, the excess tends to be 
absorbed in the thermal bath if the particles thermalise with it. However, in 
the range of parameters in which the realignment plays an important role, the 
Primakoff interaction is rather inefficient in thermalising the pseudoscalar conden- 
sate. Moreover, the process + e ± — > 7 + e ± is exponentially suppressed at high 
temperature because the energy of the incoming electron E has to be high enough 
to produce a photon whose plasmon mass m 7 is of the order of the temperature. 
This threshold implies 2m^E > m 2 , and thus E > aT 2 /m^, which is disfavoured 
especially in the case of small m^. This threshold could be overcome in a three- 
body interaction with a photon more in the initial state. However also this process 
is suppressed at high temperature. Since the FF term is a total derivative, the 
two-photon interaction involves a derivative of the pseudoscalar field whose only 
component for the condensate is do<f) ~ m^cf). Again, for small-mass pseudoscalars 
the interaction is suppressed by powers of m^/T, because all absorption amplitudes 
for the zero mode are proportional to and thus the probability to m|. 

A more subtle effect involves the presence of primordial magnetic fields, and 
could in principle absorb the excess of pseudoscalar DM pointed out in the pre- 

— * 

vious section. We have mentioned that, in presence of a magnetic field -B C xt! a 
pseudoscalar particle and the component of the electromagnetic vector potential 

— * 

parallel to B cxt mixes, because of the two-photon coupling (1.24). The two fields 
can therefore oscillate into one another. The mixing allows a pseudoscalar sur- 
rounded by a magnetic field to have effectively the same interactions as the pho- 
ton, but with amplitude suppressed by the mixing parameter XcS which is usually 
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small. These photon-like interactions should therefore be added to our thermalisa- 
tion considerations, besides the Primakoff effect. Plasma effects complicate further 
the problem, because the mixing parameter can be substantially enhanced during 
its thermal evolution. More details about the magnetic field induced evaporations 
and the derivation of the related quantities are provided in appendix A. 

The effective mixing parameter in the primordial plasma can be approximated 

by 

A - . ; f . ■ (2-60) 

Here, B ext (T) is the modulus of the magnetic field, and D(u,T) is the photon 
damping factor, which depends on the photon wavelength oo and the temperature. 

Magnetic seeds could be generated during cosmological phase transitions [124], 
thus plausible temperatures for the birth of a primordial magnetic field could be 
Aqcd or -Sew- We call T B this temperature. Because of equipartition of energy 
- and dimensional analysis — we assume for primordial magnetic fields B ~ T 2 . 
Then, the magnetic flux conservation implies that B = Bi(Rs / ' R) 2 ~ B^T/Tb) 2 
if Bi is the initial strength of the field [125], suppressing the mixing (2.60) at low 
temperatures. 

If the damping factor dominates the mixing parameter (2.60), the magnetic 
evaporation affects the pseudoscalar condensate if 

( ?± y (ea 2 (jm_\ 1/2 T b > ! (2 61) 

VlO- 10 GeV-V \nGj \g*(T B ) J 10» GeV ~ ' 1 ' 

where B is the present value of the magnetic field. An axion cold DM con- 
densate would not be affected by the magnetic field [125], since in this case 
Qa-yy «C 10~ 10 GeV _1 and B in the nG range is very close to the upper limit 
for a primordial magnetic field [126,127]. 

In the lower temperature range, when the plasmon mass m 7 is more important 
than the damping factor, there are three different regimes for XeS- We are mostly 
interested in the resonant regime, defined by m 7 (T) ~ m^, during which the mixing 
is enhanced. In the resonant case, the magnetic field induces the evaporation of 
an ALP condensate if 

1 {9^f > 1 

2 rm^H^ 
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Figure 2.10: Region of the parameter space that would be affected by the 
resonant pseudoscalar-photon mixing, leading to the evaporation of the DM con- 
densate, in case of primordial magnetic field of present strength B. 

where r = dlog {mj 2 ) /dlog(T) and the subscript res denotes quantities evaluated 
at the time of the resonance 7j res . The region of the ALP parameter space where the 
resonant enhancement can in principle affect the primordial population is plotted 
in figure 2.10, as a function of and the product g^B, where B is the today 
value of the primordial magnetic field. The two different behaviours in the curve 
are due to the resonance happening during the two different regimes of m 7 , which 
depends on the plasma electrons being relativistic or not. 

If we consider primordial magnetic fields originating at the electro-weak phase 
transition, then in order to obtain effects on the ALP population in a region 
not constrained by HB stars we would need for the field an average intensity of 
B ~ 3 /iG at present time. Cosmological observables presently exclude effective 
values of the intergalactic magnetic field larger than few nG [126,127], which leaves 
not too much space for the resonant magnetic field evaporation, at least in the ALP 
unconstrained region. 

2.4 Pseudoscalar decay 

The freezing in of the decay makes the pseudoscalars regain thermal contact with 
the bath. The decay temperature Td can be obtained from tH = 1, so we need to 
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know the energy density of the universe in order to calculate H . If pseudoscalars 
do not dominate the universe energy budget when they decay, T d is 



T« - ( 7 9 * l ) (^Vi ^ MeV . (2.63) 

gl 1 (T d ) ^10- 7 GeV-V VMeV/ 



If instead the pseudoscalar energy density does dominate, which happens when <fi 
is very non-relativistic and far out of equilibrium, the decay temperature is instead 

^ TF, ( I ) 4 ' 3 f^7) 5/3 MeV . (2.64) 

This is typically larger than the previous case, since a matter dominated universe 
expands more slowly than a radiation dominated one. Since the universe becomes 
radiation dominated after the decay, the temperature in equation (2.63) gives the 
correct order of magnitude for the reheating temperature. 

If the decay temperature Td is smaller than 3m^,, the pseudoscalar reaches 
thermal contact with the bath being non-relativistic. From equation (2.47) and 
(2.53), we know that n eq (T d ) <C n^(T d ), and thus the Boltzmann equation (2.51) 
becomes 

n + 3Hn = . (2.65) 

T 

The solution is easily found to be n(t) = n(t d ) (R d / R) 3 exp (— (t — t d )/r), where t d 
and R d are respectively the time and the scale factor when T = T d . The number 
density rapidly decreases because of the decay, and two photons per decaying 
pseudoscalar are created. 

On the other side, if T d > 3m^, n eq (T d ) can not be neglected, and the ALP 
can regain the thermal abundance — in this case it would be given by equation 
(2.41), and so higher than n^,(T d ) — thanks to the inverse decay process, 77 — > 0. 
The rate for this process is 

r^.I^tMw, (2 . 66) 

where (m^/ou) is the thermally averaged time dilatation factor, being u > the 
energy of the outcoming pseudoscalar. Decay or inverse decay is only kinemati- 
cally allowed if < 2m 7 . For T » m e , we have m 7 ~ T while for T <ti ra, 
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Figure 2.11: Axion decoupling and recoupling (C~ = 1.9, in the light yellow 
region it is also C e = 1/6). Thick solid line: freeze-out of Primakoff process. 
Medium solid line: coupling and freeze-out of the Compton process. Thin solid 
line: recoupling of inverse decay. In the yellow shaded region, axions are in thermal 
equilibrium. The light yellow region is relevant only if the Compton process is 
effective. The dashed line denotes m a = 3T, on the left of it axions are non- 
relativistic. The vertical lines delimit the e + e~ annihilation epoch. 



we have m 7 T, thus the decay/inverse decay channels open up not far from 
~ max{T, m e }. In this case, some photons are subtracted from the ther- 
mal bath to make the pseudoscalars regain the equilibrium distribution. The 
pseudoscalar population then follows the equilibrium distribution, and becomes 
Boltzmann suppressed, basically disappearing from the bath, when the tempera- 
ture drops below m^. This equilibrium decay of relativistic pseudoscalars, together 
with the previous case of non- equilibrium decay of non- relativistic particles, should 
be well kept in mind as they will be a key point to understand the topics of chap- 
ter 4. 

Focusing on the axion case, a direct coupling with the electron can be relevant 
or not according to the model. Since Tp, r and eventually the Compton interaction 
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Figure 2.12: Axion number density n a after e + e~ annihilation from numerically 
solving the Boltzmann equation until T = m e /10. The equilibrium density n e ^ 
is defined in terms of the photon temperature. Solid line: only Primakoff process 
(C 7 = 1.9). Dashed line: Primakoff and Compton (C 7 = 1.9, C e = 1/6). 

rate Tq ~ aC^n e / f a 2 are unambiguously determined as a function of f a or m a 
once a model is chosen, it is possible to plot the freezing out and freezing in 
temperatures, like in figure 2.11, where the parameters chosen are C 7 = 1.9 and 
C e = 1/6. The yellow area represents the range in which the axion interactions 
are efficient, while below the dashed line the axion is relativistic. Therefore, the 
inverse decay is active in the yellow area delimited by the dashed line and the 
thin solid one, along which T 77 ^ a = H. From this figure it is clear that the 
lower is f a , the higher m a , and thus the later the axion decouples, in particular 
above 10-20 keV this never happens. In figure 2.12 we plot the ratio between 
the axion number density and its equilibrium value, both at T = m e /10. It is 
obtained solving numerically the set of Boltzmann equations for axions, photons 
and electrons. Including the Compton interaction, the resulting yield is higher, 
for the axion, having an additional interaction channel, decouples later so it its 
abundance is less diluted by the factor (2.21). The range mostly affected by the 
further coupling is m a = 0.5-10 eV. We can see in figure 2.11 that is where lies 
the Compton light yellow region, and from this plot we can better appreciate 
the decrease in the decoupling temperature. We will deal again with figures 2.11 
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and 2.12 in chapter 4, where the Boltzmann equations will be discussed. 

We have shown in this chapter how a population of pseudoscalars can arise and 
then disappear in the early universe. Now it is time to develop its consequences. 



CHAPTER 3 



Signals from the sky: relic decay photons 

Each pseudoscalar decay produces a couple of photons, as we have seen in the 
previous chapters. The fate of these photons depends mainly at which stage of the 
universe evolution the decay takes place. 

The many charged particles in the primordial plasma efficiently scatter the 
photons propagating among them. The universe is therefore optically opaque until 
it becomes cold enough to permit the electrons which survived the e ± annihilation 
to combine with nuclei to form neutral atoms. This event, known as recombination, 
happens rather late in the history of the universe. The temperature has to cool 
down to a value around a couple of orders of magnitude below the binding energy 
of hydrogen and helium atoms. This delay trend is very peculiar of combination 
events in cosmology, as it affects also primordial nucleosynthesis. The reason is 
the overwhelming number of photons compared with the number of baryons. The 
baryon-to-photon ratio is measured to be [86] 

77= ^ = (6.23 ±0.17) x 1(T 10 . (3.1) 

n 7 

The number of photons whose energy is above the photo-dissociation threshold 
therefore can be easily larger than n^. The combination event has to wait behind, 
until this high energy tail of the photon distribution empties because of the cooling. 

The recombination is a crucial event for our understanding of the history of 
the universe. Just after this epoch the cosmic microwave background (CMB) is 
released. Most of the information we have about the early universe is extrapolated 
from what was imprinted on the CMB around the recombination time. 
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In this chapter we address the effects of the decay photons injected after the 
universe becomes transparent or immediately before. In the first case they can 
freely propagate and in principle be detected, while in the other they should leave 
an imprint on the CMB. We will analyse also the effect of the pseudoscalar decay 
on the ionization history of the universe. 

3.1 Spectral distortions of the cosmic microwave 
background 

The photon distribution is maintained in thermal equilibrium because of Compton 
scattering, double Compton scattering and bremsstrahlung interactions with the 
residual electrons. This means that if a perturbing event happens, the photon 
distribution can evolve again to the equilibrium condition, but only if these inter- 
actions have enough time to operate. Compton scattering 7 + e~ -B- 7 + e~ has 
the fastest rate among these processes, T cc a 2 , but it can only lead to kinetic 
equilibrium. Since it conserves the photon number, it can not erase a chemical 
potential in the distribution. Double Compton scattering 7 + e~ B- 7 + 7 + e~ and 
bremsstrahlung e~ + p + -B- + p + + 7 are comparatively slower, r oc a 3 . They 
change the number of photons, and thus they permit to achieve the full thermal 
equilibrium condition, i.e. a Planckian spectrum [128]. 

FIRAS measured the CMB spectrum in the range 2-21 cm -1 , and found it to 
be is a perfect black-body within very small experimental errors, only an (9(1CT 5 ) 
deviation is allowed by data [129]. Thus, any perturbation of the photon equi- 
librium condition must have been smaller than the measurement uncertainties, 
or the perturbing event must occur long enough before the bremsstrahlung or 
double-Compton cease to be effective. 

Because of the very small value of rj, the bremsstrahlung process is subdom- 
inant with respect to double Compton until e + e~ annihilation. The exact value 
of the freezing out temperature for each photon-electron interaction depends on 
the wavelength. Double Compton scattering is effective at all wavelengths until 
the temperature drops below T DC ~ 750 eV at £ DC . At temperatures lower than 
Tdc double Compton and bremsstrahlung can still produce or absorb photons, 
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but only if their frequency is u < T. Compton scattering redistributes pho- 
tons along the whole spectrum and decouples later, at Tq ~ 25 eV. Therefore, if 
a photon injection is not completely reabsorbed before t DC , the resulting CMB 
spectrum follows a Bose-Einstein distribution at high frequency and a Planckian 
one in the frequency range where double Compton and bremsstrahlung are still 
active [128]. The residual degeneracy parameter is constrained by FIRAS to be 
< 0.9 x 10' 4 [129]. 

After tnc, the photons follow a random walk path scattered by the electrons. 
Since the average photon energy is very small compared with the electron mass, 
the electron recoil is negligible and photons simply bounce off in random direc- 
tions. Photons emitted at this stage distort the overall photon spectrum in a 
very peculiar way. Electrons rapidly thermalise with the non-thermal population 
of photons and their effective temperature increases. However, as we said, CMB 
photons cannot gain energy efficiently out of them. The CMB spectrum is however 
slightly influenced by the heated electrons. They dissipate the exceeding thermal 
energy pushing few photons towards higher frequency and this imprints a typical 
pattern on the CMB spectrum. This phenomenon is described by the Kompaneets 
equation [130, 131] 



<9n 7 1 d 
dy x 2 dx 



i / & n -t . . 2 



x ( -7^- + n 7 + n 7 



(3.2) 



The variable x = oo/T e is the ratio between the photon frequency u and the electron 
temperature T e . The parameter y is given by 

dy = Te ®~ T ® a T n e (t)dt , (3.3) 

where n e is the electron number density, T the photon temperature and <7t is the 
Thomson cross section. A stationary solution of equation (3.2) is 

because of the number conserving Compton interaction. If T e 3> T, the Kompa- 
neets equation simplifies [130], 

dn -r = 1 / 35) 
dy x 2 dx dx 
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and if the perturbation to the equilibrium is not too big the approximated solution 
gives 

A™ xe x ( e x + l \ 



n 7 



whose limiting case are — 2y for i < 1 and yx 2 for x ^> 1 [130, 131]. Computing 
the exact spectrum requires to calculate the effect of the decay on T e and nume- 
rically integrate equation (3.2). Like in the /i case, the observational constraint 
is very limiting, \y\ < 1.5 x 10~ 5 [129]. Of course, all of this depends on the 
rate of photon-electron interactions. At temperature around Tdec = 0.26 eV, most 
electrons and protons are combined into neutral hydrogen, the universe becomes 
effectively transparent and the CMB decouples. 

The photon injection due to pseudoscalar decay is an example of a perturbation 
strictly constrained by FIRAS data. The non-thermal contribution to the photon 
spectrum from decay is a peak centred at energy u ~ m^/2. Assuming a thermal 
population, from equation (2.54) we have 

PjL > 0.007^ , (3.7) 

where we have taken g^s = 3.9. The energy injection is always greater than 10 _5 p 7 , 
since the ratio m^/T > 1, and thus it is potentially in radical conflict with ob- 
servations. To conventionally escape this limit only a small fraction of the energy 
injection should influence the photon distribution, thus the decay has to be just 
started when CMB is released or to occur long before recombination. However, 
we will see in the following that there are other cosmological data sets which ex- 
clude pseudoscalars able to clear the CMB distortion hurdle through late or early 
enough decay. These further limits largely overlap with the CMB distortion bound. 
Therefore, here we consider sufficient only a rough estimate of the pseudoscalar 
parameter space which is excluded by CMB spectral distortions. A cosmological 
scenario with g*(T fo ) > O(10 3 ) would permit to evade unconventionally this CMB 
constraint. Also the dilution of the pseudoscalar density due to some consistent 
entropy injection during the epoch in which is decoupled from the thermal bath 
can in principle relax this bound. However, this has to happen much before pri- 
mordial nucleosynthesis, otherwise it can be constrained by the same arguments 
which will be applied to pseudoscalars in the next chapter. 
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Figure 3.1: Left: Photon /i parameter after axion decay (5 = 1). The obser- 
vational upper bound is indicated. Right: The contours of log 10 \/J.\ in the 5-m a 
plane where black/red corresponds to positive/negative values of /x. The thick 
black line corresponds to the boundary /i = 0. The shaded area is excluded. 

In figure 3.3 — which gives a summary of this chapter — the light green regions 
labelled CMB /j, and CMB y corresponds to Tbc > T&>Tq and Tq > Id > ^dec 
respectively. The influence of later decays requires a different treatment and it is 
discussed in the following sections. We always assumed that the amount of photons 
injected is always too large to be thermalised, since we are not considering here the 
unconventional scenario of pseudoscalar dilution just discussed. Note that these 
bounds are somewhat conservative. The decay is not an instantaneous event, so a 
significant amount of energy can be released after T&. This is especially true for 

> keV, for in this case the pseudoscalar energy density dominates the universe 
before the decay and the subsequent perturbation of the photon spectrum is huge. 

In the axion case, we have numerically solved the differential equation for the 
evolution of the degeneracy parameter \x of the photon distribution [1]. Assuming 
\x small, its evolution is governed by the equation [128], 

|^_, ( r DC + r B ), (3.8, 

where Trjc and I~b are the double Compton and bremsstrahlung rates for decreasing 
/i, which are the inverse of the relaxation times for the two processes. The rate- 
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of-change due to the axion decay photon injection is given by 

dfi a = 2_ /j^^ _ J^cin„\ 

dt 2.14 Vp 7 dt n 7 dt J ' 
where dn a /dt = —n a /T and dp a /dt = m a dn a /dt. Our results are shown in fi- 
gure 3.1, where we plot the final /i value as a function of m a and 5 = C 7 /1.9 for 
hadronic axions. In the left panel we fixed C 1 = 1.9, the value in the simplest 
KSVZ model, and in this case we found 

m a > 8.7 keV at 95% C.L. . (3.10) 

This is a robust bound, since /j, is a steep function of m a . The CMB distortion 
effect depends sensitively on the axion-photon interaction strength for 5 < 1 (right 
panel of figure 3.1). Generally the spectral distortions get larger for smaller C 1 
at a given m a , because if the decay happens later the photon distribution is less 
protected against distortions. For large C 7 , the final \i changes sign from negative 
to positive with increasing m a . For 5 < 0.1, fi is always positive since axions 
decay non-relativistically, thus the energy injection is more important than photon 
number, see equation (3.9). Of course, because of the sign change in /i some fine- 
tuned cases exist where the final \i can be accidentally zero. 



3.2 Ionization history of the universe 

CMB is released at around T dec , whose corresponding redshift is Zd ec ~ 1100, 
when most of the recombination has already taken place. Before recombination, 
the number of ultraviolet photons was too high to have long-lived atomic bonds, 
which were broken just after their formation. 

Around z TCC ~ 1500, when the universe temperature was T ~ 0.35 eV, helium 
and hydrogen nuclei start to capture free electrons to form neutral atoms. Helium 
recombines before hydrogen because of its higher atomic number and thus larger 
Coulomb potential. The universe is still opaque after helium recombination, since 
He constitutes just around one fourth of the baryonic density, and the CMB has 
to wait the hydrogen recombination to break free [132]. 

To exactly compute the ionization history of the universe it is necessary to 
solve a set of coupled Boltzmann equations for photons, free electrons and atoms 



3.2 Ionization history of the universe 



63 



in their different energy levels. The result of this calculation is usually given in 
terms of the fraction of free electrons, or ionization fraction Xi OQ , as a function 
of redshift. The recombination of the hydrogen freezes out at z ~ 800 because 
of the expansion of the universe, leaving a residual ionization fraction of order 
Xion ~ O(10 -4 ). The universe fully ionizes again much later, between redshifts 6 
and 10, presumably due to ultraviolet emission from the first galaxies. The details 
of the reionization process are still not well understood [133]. 

The slight imprint that the free electrons leave in the CMB through Thomson 
scattering, for instance in the polarization, gives us information about the history 
of recombination and reionization. The optical depth for CMB photons is one of 
the parameters that can be measured from the CMB multipole analysis and it is 
defined to be 



where all the quantities in the integral are expressed as a function of redshift. 

The WMAP7 measured r opt after decoupling to be 0.088 ±0.015 [86]. A factor 
0.04-0.05 of this can be attributed to a fully ionized universe up to about redshift 
6, which is supported by the absence of Ly-a features in quasar spectra. The 
origin of the remaining fraction, r 6 = r opt (6, z dcc ) ~ 0.04, is still uncertain and 
leaves some space for ALP decay. 

The photons produced by ALP decay after decoupling are free to propagate 
since there are almost no free electrons to interact with. However, ultraviolet ra- 
diation can interact with atoms and photoionize them. The universe is indeed very 
opaque to ultraviolet radiation. This prevents us from detecting decay photons 
in the 13.6-300 eV range, but it does not mean that we can not constrain them, 
since the photoionisation triggers an increasing of the ionization fraction of the 
universe. This argument does not hold for axions because in this mass range they 
decay before recombination. If we assume that each decay photon ionizes only one 
H atom immediately after the ALP decay, which is a first rough approximation, 
the number of ionizations per unit time can be estimated to be [134] 




(3.11) 
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Figure 3.2: The effect of ALP decay on the ionization fraction, plotted in 
function of the reshift z. The standard case with no decay is plotted with the 
solid line. Left: decay of a 100 eV mass ALP with r = 10 18 s, 10 20 s, 10 22 s, 10 24 
s and 10 26 s. Right: decay of a r = 10 22 s ALP with = 1000 eV, 500 eV, 100 
eV and 50 eV. 

If we now multiply this quantity by a typical time scale [134] 

t H = 1/H(z) ~ 2.4 Myr [501/(1 + z)f /2 , (3.13) 

where we have neglected fl r and fl^ in H(z), we get a conservative estimate of the 
degree of ionization induced by the ALP decays up to a certain redshift. ALPs with 
100 eV mass and ~ 10~ 13 GeV -1 would have produced an ionization comparable 
with the standard residual value 10 -4 already at high redshifts, z ~ 500. Since 
the ALP-induced ionization grows in time as (1 + z)~ 3 ^ 2 , it shows a potentially 
interesting effect. We could even think that ALPs close to these parameters may 
be responsible for the full reionisation of the universe. However, the (1 + z)~ 3 / 2 
dependence is too soft — reionization seems to be a much more abrupt process, 
usually parametrised to be almost a step function in x ion (z) — and we can only 
check if ALPs provide a r 6 compatible with observations. 

In order to obtain a more detailed constraint, we have calculated the ionization 
history of the universe in the decaying ALP cosmology [2,3] by introducing the ALP 
ionizations in the recombination code RECFAST [135]. In figure 3.2, we plot the 
ionization history of the universe in function of z for several examples of ALP mass 
and lifetime. We set the ALP fraction to to account for all the DM, O dm /j 2 = 0.11, 
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and made the conservative assumption that each photon emitted during the decay 
can ionize only one atom. On the left panel the mass is fixed, = 100 eV, and 
the lifetime varies. Of course, for longer lifetimes the ALP ionising effect appears 
later. On the right panel the lifetime is instead fixed to 10 22 s, and the mass varies 
from 50 eV to 1 keV. The higher the mass, and consequently the energy u of the 
emitted photons, the less efficient is the ionising effect. The one-electron atom 
photoionisation cross section is suppressed for very high energy photons [136], 



where Z is the atomic number, E\ s = 13.6Z 2 eV the energy of the Is state, 
do = (am e )"' = 5.292 x 10~ 9 cm is the Bohr radius and u the photon energy. 

To scan the ALP parameter space, we computed the optical depth in the inter- 
val z = 6-100, requiring it to not exceed r 6 . We made two different calculations, 
assuming the ALP thermal abundance in [2], and secondly that ALPs constitute 
the whole DM in [3] . Our results are excluding the light green region labelled x- loIi 
in figure 3.3, where the thermal origin of ALPs is considered. A similar result 
was obtained in [3]. This bound would increase up to one order of magnitude at 
the largest masses for which ionization is effective, < 300 eV, if we assume 
optimistically that all the energy of the emitted photons can be converted into 
ionization. The ionization history constrains ALP lifetimes much longer than the 
age of the universe, r > 10 24 s, which means that only less than one ALP out of 
ten millions can decay. The effect of the decay of a large population of particles 
has catastrophic effects, but only extremely small perturbations to the standard 
cosmological scenario are allowed. 



We have seen in the previous section that ultraviolet light is absorbed by atoms. 
If the photons from ALP decay do not lie in the UV region described above, the 
universe is transparent to them and they can in principle be detected. The decay 
photons should be manifest in spectral observations as a peak at energy u> = 
broadened by redshift. In the parameter space we can constrain, the decay happens 




(3.14) 
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at rest in the comoving frame. The spectral flux of photons produced in the decay 
of a diffuse pseudoscalar population is [31,137] 

dFE - 1 1 - (3.15a) 



dEdQ. 2nH(z) (1 + zf 

~ 2tttH \m+/2) 6XP 
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3/2' 



r \m (j> /2j 



(3.15b) 



where the subscript means quantities at present time, n^ is the putative number 
density if would be stable, and E is the energy at which the photon would be 
seen today. The photon initial energy is m^/2 and thus the redshift of the decay 
is given by 1 + z = (m<f ) /2)/ Eq. For simplicity we assumed matter domination 
neglecting 0(1) corrections due to the cosmological constant. To take into account 
the opacity of the universe to the ultraviolet radiation, the flux has to be corrected 
multiplying it by the survival probability 

P( z ) = e ~ K ^ E) (3.16) 
f z dz' 

k{z,E) = + ^ rcg(*V P h-ion(£) , (3-17) 

where E = E (l+z), Uh is the hydrogen number density and <7 p h-ion is the hydrogen 
photoelectric cross-section given by equation (3.14), with E u = 13.6 eV and Z — 1. 
The contribution from helium photoionization is not relevant. Once we corrected 
the spectrum by this factor, we compared it with the extragalactic background light 
(EBL) spectrum reviewed in [40]. We corrected also the microwave range of the 
spectrum subtracting the CMB contribution, which is known extremely well after 
the measurement of COBE-FIRAS. The EBL spectrum does not present sharp 
features that eventually could be recognised to be the contribution of a decaying 
particle. Therefore we require the decay photon spectrum to be lower than the 
EBL. The region of parameter space excluded by this comparison is plotted in 
dull green in figure 3.3, where it is labelled EBL. The same approach was already 
used in references [39,40] to look for axions in the optical EBL. In the part of 
parameter space that can be constrained in this way, the ALP abundance would 
be the minimum considered n^/n 7 = (3.9/106.75)/2 ~ 0.019. Since the thermal 
ALPs would provide too much DM for > 154 eV, in our analysis we found more 
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convenient to assume that ALPs provide the right amount of DM above this mass, 
while below it the standard thermal abundance is considered. In this way, we 
avoid to duplicate the limits in the region which is already excluded, and provide 
a bound that only implies ALPs to constitute the DM. Of course this requires 
some non-standard dilution of the ALP number density by additional degrees-of- 
freedom above the electroweak scale. The bound gets severely degraded in the 
= 13.6-300 eV range, where not only absorption is very strong but also the 
experimental data are extremely challenging and the EBL spectrum has only an 
upper bound estimate. Of course these facts are closely related. 

In order to gain sensitivity, a more effective strategy to search for decay photons 
is to examine the light emitted in galaxies and large scale structures. There, the 
dark matter density is above the average, providing therefore an enhanced signal 
for the decay products. However, also the backgrounds are above the EBL, but 
anyway the decay photons should stand above the background appearing like a 
peak. Searches of axions in the visible have been presented in [38, 39] and more 
recently in [41]. The limits from these references are plotted in figure 1.2 in the 
blue band labelled Telescope. In [138, 139] sterile neutrino dark matter decays in 
a photon plus an active neutrino were searched in X-rays with the same methods. 

We followed these references, rescaling the results for the ALP decay mode when 
exploiting neutrino searches, and abundance when the thermal ALP yield gives 
Q$ < f^DM- The exclusion bounds we obtained are plotted in figure 3.3, labelled 
respectively Optical and X-Rays. We have also plotted in figure 3.3 the limit 
obtained by galactic line searches in the 7-ray range for decaying DM, labelling it 
7-rays. We took the data from [140] and again we rescaled them for the ALP decay 
mode. Anyway, considering that the region constrained by 7-ray observations 
is well above the Tf D > mpi line, it makes probably not too much sense in our 
thermal-ALP picture, but we showed it anyway for the sake of completeness. In 
this 7-ray region, a very suggestive yet tentative claim of a detection at uo ~ 130 
GeV appeared recently [141]. Unfortunately, we can not relate it to decaying 
pseudoscalars, since it can only be attributed to annihilating DM. 

In these spectral analyses, the monochromatic emission line caused by the 
decay — or annihilation — is reconstructed, and its intensity tracks the density 
profile of the galaxy or the cluster taken into account, which can be measured 



68 



3. Signals from the sky 






W~ 

% 


— i, — ^ ^ — 

\ ', 
v 

HB \ ^ 


CMB y 


\\ . 


- 

1 1 1 1 


, i , i ) 


excluded ' \ V 1 

**-\ \\ 


CMB n 
1,1, 


1,1,- 



4 6 
Log 10 [eV] 



10 



Figure 3.3: Summarising picture for this chapter in the m^-T parameter space. 
All the bounds are described in the text. In light green are plotted the bounds 
related to CMB distortion, respectively labelled CMB y, CMB /x and x- lon . In 
dull green are the limits from the lack of direct observation of decay photons. 
They are EBL and Optical, X-Ray and 7-Ray. We considered the ALP 
abundance provided by the thermal production mechanism for all the bounds 
in all the parameter space but direct decay photon detection in the > 154 
eV range. Here, ALPs would be overproduced with respect to the measured 
DM abundance, therefore SIdm is taken for ALP abundance in the direct photon 
detection limits. The orange region labelled DM covers the exclusion bound 
provided by ALP overproduction. 
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through lensing observations. The high resolution of modern instruments permits 
to extract different spectra from different zones of an object. The emission line 
are thus searched privileging the higher density regions. 

The radiation produced by ALP decay being proportional to the relic ALP 
density would be linearly sensitive to the particle content above the EW scale 
through <7*(Tf ). The constraint on relaxes as a/T+ new dof / 106.75 with the 
new thermal degrees of freedom with masses between the EW and T fo in the 
< 154 eV range. The new degrees of freedom do not affect the direct decay 
photon detection bound for > 154 eV, since in this mass range we assumed 



CHAPTER 4 



Diluting neutrinos and nucleons: 
limits from dark radiation and big-bang 

nucleosynthesis 



If the pseudoscalar decay happens long before matter-radiation decoupling, the pri- 
mordial plasma has enough time to regain thermodynamic equilibrium. This does 
not mean that is not possible to observe any consequence of this event. The decay 
modifies the photon abundance, and leaves unaltered the baryon and the neutrino 
abundances — if neutrinos have already decoupled from the bath. The baryon-to- 
photon ratio rj and the effective number of neutrinos N e g oc p v j ' p 1 are effectively 
diluted, and they can be measured through the analysis of CMB anisotropics and 
the observation of primordial elements abundance. 

A fast estimate of the entropy increment due to the decay is a useful guideline 
to understand the topics of this chapter. Since neutrinos are not affected by the 
decay if already decoupled, it proves convenient to normalise abundances in terms 
of the neutrino temperature which redshifts as T v oc RT 1 . So we define 



The initial condition we use is T 7j0 = T ufi = T^ = 2 MeV, when neutrinos have 
decoupled and e + e~ annihilation has not yet begun. A first stage is defined by 
T 7; i ~ m e /10, when electron-positron annihilation is over. Entropy conservation 




and 




(4.1) 
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during e + e annihilation implies 

73 Q + 2 + =Ti(2 + &), (4.2) 

where 7/2 and 2 are the e + e~ and 7 entropy degrees of freedom. Because of the 
initial condition we have chosen, 7o and £ are equal to 1. Considering the axion 
case, from figure 2.11 we see that we are describing an epoch far below the dashed 
line, which means that they are relativistic in the mass range we are interested. We 
assume kinetic equilibrium, even if axions decouple during the e + e~ annihilation. 

If axions or ALPs decouple before e + e~ annihilation, like neutrinos, they are 
not affected by the entropy released in the annihilation. Photons are heated by 
the standard amount, thus 7i = 11/4, while the product £T is conserved. On the 
other hand, if they decouple during or after e + e~ annihilation they are heated, 
sharing some of the electron entropy, and thus we have 

71 = 2KTT2) • (43) 

If pseudoscalars are fully coupled during e + e~ annihilation, then £1 = 1 and 7i = 
13/6. In general we have 11/4 < 71 < 13/6. Figure 2.12 plots the value of £ for 
the axion at 7\. 

The eventual decay of a population of pseudoscalars makes the entropy they 
carry to be transferred to the photon bath. In the case that never leaves ther- 
mal equilibrium, we easily find that the final photon abundance is 13/11 times 
the standard value. In the axion case, we have pointed out that this happens for 
m a > 20 keV. In general the entropy transfer to photons depends on two parame- 
ters, the initial pseudoscalar abundance, parametrised by £, and the effectiveness 
of inverse decay when axions and ALPs become non-relativistic. 

If the decay freezes-in when T 7 > m^/3, the pseudoscalar temperature catches 
up with the photon one which has to decrease because of energy conservation. In 
this process, radiation is simply shuffled from one form to another and comoving 
energy remains conserved. Using again neutrinos as a ruler, in our units the 
energy density is p 1+< p/p u oc (2 + ^ 3 )Ti^ 3 before the recoupling. After recoupling 
at T*, the pseudoscalar temperature catch up with the photon one, thus we have 
p 1+( p/p u oc (2 + l)7^ 3 , and, because of energy conservation, the photon abundance 
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is reduced to 




3/4 



(4.4) 



Comoving entropy increases by the factor [(2 + l)/(2 + £i)]%/7i- At the later 
temperature T 2 , when pseudoscalars become nonrelativistic and their abundance 
gets Boltzmann suppressed, the decay transfers the pseudoscalar entropy to the 
photon bath, heating it according to 



Putting all together, the final photon heating by pseudoscalars recoupling and 
adiabatic decay is 



which is also the ratio of the entropy density after and before the decay. We refer 
to this case as equilibrium decay, since it happens in a local thermal equilibrium 
condition. 

If pseudoscalars recouple non-relativistically and they dominate the energy 
budget of the universe just before recoupling, their out- of- equilibrium decay has 
more dramatic consequences. The out-of-equilibrium decay can produce a vast 
increase in the entropy density, because the energy which is first stored in the 
mass of a non-relativistic particle is suddenly released as radiation. This case is 
illustrated by an analytic approximation to the entropy generation [76] 



where (gj s ) denotes an average over the decay time and = n^/s. This formula 
is valid if the decay produces relativistic particles, and the energy densities of 
all the species other than are negligible before the decay. The energy density 
scales as p oc R~ 3 for non-relativistic species, while p oc i?" 4 for radiation. A 
pseudoscalar dominated universe requires non-relativistic particles which are long- 
lived enough to survive until the pseudoscalar-radiation equality. This condition is 
more easily met in the intermediate mass range for a given coupling constant: very 



2T 2 = (2 + 1)T* . 



(4.5) 




(4.6) 




(4.7) 
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Figure 4.1: Photon density increase in our modified cosmology as expressed by 
S2/S1 for C 7 = 1.9. Solid and dashed lines stand for hadronic and non-hadronic 
axions respectively In the latter case C e = 1/6. The thin blue lines show the 
value if we assume that no entropy is generated in axion decay. The vertical line 
signals the end of e + e~ annihilation epoch. 

heavy pseudoscalars decay early because of the ml 3 dependence of the lifetime, 
while very light particles becomes not-relativistic late. 

In figure 4.1 we show using black lines the resulting S2/S1 from a numerical 
solution of the set of Boltzmann equations for the axion phase space distribution 
and the neutrino, photon and electron abundances as a function of m a in the case 
0^ = 1.9, with and without the electron coupling. The blue thin lines represent 
just the increase in the entropy density due to reshuffling of entropy between 
axions and photons, which follows equation (4.6), ignoring the entropy generated 
in the out-of-equilibrium decay. It is interesting to observe how dramatic is the 
effect of the axion domination on the comoving entropy density. Starting from the 
very right, in the high mass range, we can observe how negligible is the effect of 
the decay, since in this case the axion population has already disappeared from 
the bath when T = T\. The value of S2/S1 increases moving to the left, as the 
equilibrium decay happens more and more after the electron-positron annihilation. 
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The maximum value allowed in this situation is provided by equation (4.6) and 
is 3/2. Continuing to the left, especially if no direct coupling with the electron 
is present, axions experience a period of decoupling followed by recoupling due to 
inverse decay. The value of S2/S1 decreases because £1 is increasingly less than one 
in equation (4.6). Finally, the out-of-equilibrium range begins, and the black and 
blue lines diverge. Here, ^ tends to zero, and indeed the blue lines approach two. 
But the actual value of S 2 /S 1 is affected by entropy generation, and it is highly 
boosted if the decay happens late enough for the axion to dominate the universe. 
Again, the additional coupling to electrons makes the transition from equilibrium 
to non-equilibrium decay shift towards lower mass. 



From the previous discussion, we have seen that the pseudoscalar decay can in- 
crease T, i.e. the abundance of photons relative to neutrinos. However, the cos- 
mological point of view is a bit different: since all cosmic parameters are defined 
relative to the observed CMB properties, this means that effectively the neutrino 
abundance is reduced. Precision observables measure the cosmic radiation density 
at decoupling, and this neutrino dilution should make a difference. 

To extract cosmological information on the radiation density at CMB decou- 
pling, in [1] we analysed the usual 8-parameter standard ACDM model described in 
reference [122], extended in two ways. We allowed the effective number of neutrino 
degrees of freedom to vary, assuming a flat prior on the interval < N cS < 3.0. 
The radiation energy density is traditionally expressed as 



where N eS is the effective number of thermally excited neutrino degrees of freedom. 
The standard value is N cS = 3.046 instead of 3 because of residual neutrino heating 
by e + e _ annihilation [142] but, given the experimental uncertainty, we neglected 
this tiny correction. In addition, we allowed for neutrino masses, assuming a 
common value m v for all flavours and a flat prior on < Q u /Q m < 1. The other 
parameters and their priors are identical with those provided in reference [122] and 
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Parameter 



Prior 




0.4-1.0 



0.005-0.1 



0.01-0.99 



Topt 



0.01-0.8 



ln(10 10 A s ) 



2.7-4.0 



0.5-1.5 



Table 4.1: Priors for the cosmological fit parameters [122]. All priors are uniform 
in the given intervals. 

are listed in table 4.1. They are the cold dark matter density f2 cdm , the baryon 
density fib, the present day normalised Hubble parameter h, the optical depth to 
reionization r opt , the amplitude of the primordial scalar power spectrum A s , and 
the scalar spectral index of the primordial fluctuation spectrum n s . Moreover, 
we use the same set of cosmological data, which are the WMAP 7-year CMB 
measurements, the 7th data release of the Sloan Digital Sky Survey, and the Hubble 
constant from Hubble Space Telescope observations. 

Marginalizing over all parameters but m v and N c q we find the 2D credible 
regions shown in figure 4.2. Marginalizing in addition over m v we find the limits 



These limits are very restrictive, because on present evidence cosmology actu- 
ally prefers extra radiation beyond N e g = 3 [86,143-146]. The favourite value is 
N cS = 4.34l°;j| considering only WMAP7 data, and N cS = 4.78lJ;fg if also SDSS 
observations are taken into account [144]. 

The PLANCK satellite, currently taking CMB data, is expected to boost pre- 
cision determinations of cosmological parameters. It will measure the cosmic ra- 
diation content at CMB decoupling with a precision of about AN eS = ±0.26 or 
better and thus will clearly decide if there is extra radiation in the universe [147]. 
If it finds convincing evidence for extra radiation, much in cosmology will have to 
be reconsidered besides our axion-ALP limits. 





2.11 at 99% C.L. 



at 95% C.L. 



(4.9) 
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4.1.1 ALP bounds 

At the beginning of this chapter we have described some limiting cases to un- 
derstand the matters related to dilution and photon injection. The key quantity 
was T . But considering the large ALP parameter space and all the intermediate 
cases it includes, it is unavoidable to numerically solve the problem. The active 
species to consider are photons, electrons and neutrinos. Electrons acts as messen- 
gers between photons and neutrinos. Photons thermalise very fast with electrons, 
which, through weak reactions, heat the neutrinos to keep track of the temperature 
changes in the electromagnetic bath. The most relevant of these energy redistribu- 
tion processes is e + e~ — > vv. We neglect scattering processes such as e ± u — > e ± i > , 
since they cannot change the neutrino number and are less effective in transferring 
energy to the neutrino bath 1 . The energy flow from ALPs to neutrinos can then 
be modelled by a set of Boltzmann equations for comoving energies, defined as 

X, = p t R\ (4.10) 

1 Assume that electrons have temperature larger than neutrinos like in our case. Then the 
speed of the energy transfer per unit volume in scattering processes is proportional to (SE)T^T^, 
with (SE) ~ (T e — T„) a thermal-averaged energy transfer per scattering. For annihilations it is 
oc T e (T^ — T e 8 ), which is much more sensitive to nonequilibrium situations. 
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and the ALP phase space distribution function f{K^) as a function of the comoving 
momentum, = Rk^, 

±f( K4> ) = -(c, + c P )(f-n, (4.11a) 
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(4.11b) 
(4.11c) 
(4.11d) 

(4.11e) 



where = \J k^ + m^ 2 is the ALP energy and = u^R. The set of equa- 
tions (4.11) describes the evolution of the comoving energy density stored in 7 
together with e ± , the three species of z/s, the ALPs and the cosmic scale factor R. 
The solution of equation (4.11a) is the distribution function / for the pseudoscalar 
particle. The collision terms for the decay and Primakoff processes are [1] 
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Since the energy stored in electrons and positrons is transferred directly to the pho- 
ton bath, we kept X 7 and X e together in equation (4.11b). The term 5p e = p e —p e /3 
accounts of the comoving energy density gain as electrons become increasingly 
non-relativistic, experiencing the transition from p e oc R~ 4 to p e oc R~ 3 . Here, p e 
and p e are the pressure and energy density of e ± . The second term in the right 
hand side reckons up the photon- ALP energy exchange, while in the last last line 
of equation (4.11b) are the neutrino-electron interactions. The equations (4.11c) 
and (4.11c) provide the evolution of neutrino comoving energies. Neutrino flavour 
influence the electron-neutrino energy exchange rate. We have separated the three 
neutrino species in two unknowns, according to the presence of not of charged 
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current interactions with electrons. The first is the case of X„ e for electron neu- 
trinos, while in the second one we have X v T for muon and tau flavours together. 
The energy exchange rate is given by the factor T eu = G 2 F T 7 , with Gp the Fermi 
constant, multiplied by c e ~ 0.68 for the electron flavour and c MT ~ 0.15 for the 
muon and tau ones, which follow from the appropriate thermally averaged cross 
section. We assume that neutrinos always have a thermal distribution, determined 
only by an effective temperature, which should be a reasonable first approximation 
and accurate enough for our purposes. We neglect the energy reshuffling between 
different neutrino species, which does not influence the total neutrino density at 
leading order. Finally, through equation (4. lie) we calculate the expansion rate 
of the universe through the evolution of the cosmic scale factor R. The initial con- 
ditions are specified at T > MeV by having all species at a common temperature 
and the ALP number density given by equation (2.53). 

For values < 10~ 7 GeV -1 , the Primakoff process is decoupled in the temper- 
ature range of interest and can be neglected. If the inverse decay is also negligible, 
we can integrate the ALP phase space distribution explicitly and directly compute 
the evolution of the number density 

li ^ = (414) 

instead of using equation (4.11a). In this way we recover the exponential decay 
law Nfi oc e~ l l T . The integral in equation (4.11b) is then 

/ f$^ (C7 + CP)(/ " /CQ) " m * n * Ri l T ■ ( 415 ) 

We have scanned the ALP parameter space in [2], and we present our results for 
N e ff in the g^-r and m^-r planes in figure 4.3. We required a standard cosmology 
scenario at temperatures below the standard matter-radiation equality T mr ~ 1 eV. 
Thus we have analysed only the part of parameter space where > T mr , in order 
to not perturb the standard picture with the decay. From equation 2.63, we obtain 
a mass lower bound for the validity of our calculations 
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Log 10 [MeV] Log 10 g [GeV '] 

(a) (b) 

Figure 4.3: Contour plots of the effective number of neutrinos N c g as a function 
of the ALP mass and lifetime (left) and of the coupling parameter and lifetime 
r (right). ALP cosmologies leading to N c g < 2.11 can be safely excluded. 



which largely overlaps with the limits obtained from CMB spectral distortions 
obtained in the previous chapter. 

The sharp features of the isocontours in figure 4.3 are due to the abrupt and 
sizeable decrease of g*s during the QCD phase transition. The relic abundance 
of ALPs depends on through the freezing out temperature, and for ALPs with 
coupling around g^ ~ 1CT 7 - 5 GeV" 1 we have T fo ~ Aq C d from equation (2.58). 
In a small range of the ALP abundance considerably changes, according if 
the decoupling is late enough for ALPs to share the entropy transferred from the 
QCD degrees of freedom or not. The larger is the ALP abundance, the larger is the 
subsequent entropy release at the time of the decay, hence the peculiar shape of 
the isocontours. In figure 4.3a we observe the same features, as the g<f, dependence 
is hidden in r. 

In figure 4.4, we present some illustrative examples depicting the evolution of 
the Xs of electrons, neutrinos and ALPs as a function of the temperature. Note 
that in all of them, when ALPs become non-relativistic, the ratio X ( f ) /X 1 rises 
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because it becomes proportional to m^/T, until the age of the universe becomes 
comparable with r. 

Figure 4.4a shows a typical case of very early decay of a massive ALP, when 
neutrinos are still partially coupled to electrons. This example depicts the ALP 
behaviour in the lower right corner of figure 4.3a, and correspondingly in the lower 
left one of figure 4.3b. Here N eS = 2.6, marginally different from the standard value 
of 3. Even if the ALP energy dominates the universe and the entropy injected du- 
ring the decay is huge, the reheating temperature is large enough for neutrinos to 
almost fully recover their thermal abundance. In general, the final value of N cS is 
related not to the total entropy released, but only to the part of it injected after 
the freeze out of the neutrino-electron interactions. In this particular region of 
the parameter space, the neutrino dilution is mainly sensitive to the ALP lifetime 
and not to m<j, or individually. The outcome of a decay earlier than t ~ 10~ 2 s 
is undistinguishable from standard cosmology, since neutrinos regain completely 
their thermal abundance. Around T ~ m e , electrons and positrons become non- 
relativistic and annihilate, heating the photon bath but not the neutrinos, which 
have decoupled. The ratio X v /X 1 therefore decreases. In this period, the photon 
temperature increases with respect to the neutrino temperature by the standard 
factor (4/11) 1 / 3 due to entropy conservation. 

The picture changes considering later decays, after neutrino decoupling: pho- 
tons and electrons get all the ALP enthropy, as shown in figure 4.4b. This is what 
happens in figure 4.3a and in figure 4.3b above the thick isocontour. The neutrino 
dilution is computable in this the ratio of the final and initial comovine: 

entropies of the photon-electron bath is given by equation (4.7). The tempera- 
ture of the electromagnetic bath increases with respect to the neutrino one by a 
factor (Sf/Si) 1 / 3 , making N cf[ = 3(5 / / S^)" 4 / 3 because also the electron entropy 
ends up in photons. The neutrino energy density is therefore strongly diluted by 
the energy gain of photons plus electrons. Note that this is mainly a function of 
m^y/r oc [g^^/m^j 1 — as we can evince from equation (4.7) — which produces 
the characteristic slope of the isocontours at long ALP lifetimes — r oc m^ 2 in 
figure 4.3a and r oc g^ in figure 4.3b. For < 1CT 9 GeV -1 , this is the only 
dependence on g^, since is constant. Another example of late ALP decay, but 
with a smaller mass, is shown in figure 4.4c. In this case we observe first the e ± 
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Figure 4.4: Examples of the evolution of the comoving energy of all species 
of neutrinos (solid), ALPs (dashed) and electrons (double-dashed), as functions 
of the temperature of the universe. All energies are normalized to one thermal 
bosonic degree of freedom. 
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annihilation, which heats photons with respect to neutrinos. A sizeable neutrino 
dilution is observable after the ALP decay. 

Finally, in figure 4.4d we show an example for which inverse decays are relevant, 
which is on the left side of figure 4.3a and on the right in figure 4.3b. As the 
temperature drops, we observe a first decrease of the ALP energy density due to 
the electrons heating the photon bath. The inverse decay channel opens around 
T ~ 70 keV and helps ALPs to regain equilibrium before disappearing at T ~ 
m<£. During rethermalization, the photon energy decreases, which can be seen 
as a slight rise in X u /X 1 . Due to this mechanism, in the small-mass and short- 
lifetime region of the parameter space we have seen that entropy conservation gives 
N c s = 3(11/13) 4 / 3 ~ 2.4. If is larger than a few MeV, the decay-in-equilibrium 
happens when neutrinos are still coupled, so N eS approaches 3. The disappearance 
from the thermal bath is governed only by the mass, and the isocontours of N eS 
exactly follow the isocontours of m^. 

The excluded region of the parameter space are determined comparing these 
numbers with the limits (4.9), taking the 99% C.L. value of N cf[ > 2.11. We are 
compelled to this conservative choice because the 95% C.L. value of N cS > 2.39 is 
just below the N eS ~ 2.4 obtained from the case of decay-in-equilibrium. Given all 
the uncertainties of this calculation, we can consider this large part of parameter 
space disfavoured but not excluded. In figure 4.3, the excluded region is above the 
thick line, which is coloured in yellow in the summarising plot for this chapter, 
figure 4.10. 

4.1.2 Axion bounds 

Pertaining to neutrino dilution, axion phenomenology is simpler than that of ALP. 
The cosmic energy density in neutrinos is modified by the factor (T^ x /T* td ) A be- 
tween the axion and standard cosmology. The results of the introduction imply 
that this ratio can be expressed in terms of the quantity T2, the modified {T 1 /T v f > 
value after axions have disappeared, and 




(4.17) 
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Figure 4.5: Radiation density during the CMB epoch (C 7 = 1.9). Solid and 
dashed lines are for hadronic and non-hadronic axions (C e = 1/6). The 95% and 
99% C.L. lower limits from Eq. (4.9) are shown as horizontal lines. 

The variation of N cS with m a for C 7 = 1.9 is shown in figure 4.5. At sufficiently 
high m a , the inverse decay process keeps thermal equilibrium during the decay. In 
this case, we have an analytical expression 



For m a > 20 keV we have £i ~ 1 and N c g reaches asymptotically the minimum 
neutrino dilution 3 (11/13) 4/3 = 2.401. At much lar ger masses, m a ~ MeV, axions 
would disappear in local thermal equilibrium before neutrino decoupling, leaving 
no trace in cosmology. Therefore, at m a ~ MeV the value of N e s shown in fi- 
gure 4.5 reaches a plateau at the the standard value N e s = 3. Comparing with 
the cosmological limits of equation (4.9) we see that even the minimum neutrino 
dilution is only barely allowed at 95% C.L. and, again, disfavoured, but not credibly 
excluded. 

For m a < 20 keV, axions decay increasingly out of equilibrium, creating entropy 
and reducing the final N e g further. Comparing the calculated N c s in figure 4.5 




(4.18) 
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with the observational limits implies 



m a > 3 keV at 99% C.L. . 



(4.19) 



4.2 Influence on big-bang nucleosynthesis 



During the primordial nucleosynthesis neutrons and protons combine to form nu- 
clei. Because of the absence of stable nuclear combinations with mass number 5 
or 8, only nuclei lighter than beryllium can form. Heavier elements will appear 
only after the formation of the first stars, where the density is high enough to 
make triple a collisions possible. We can follow the development of BBN and the 
evolution of the primordial nuclear abundances in figure 4.6. 

The ratio of nucleon abundances set the initial condition for BBN. For T ^> 
1 MeV and f < 1 s, the weak reactions 



where Q = m n — m p = 1.293 MeV is the difference between the nucleon masses, 
n n is the neutron number density and n p the proton one. Relativistic electrons 
and neutrinos are more abundant than nucleons of a factor close to the baryon-to- 
photon ratio 77, and eventual deviations of n n /n p from the equilibrium value (4.22) 
are rapidly erased at this stage. Around T ~ 1 MeV and t ~ 1 s, the weak 
reactions (4.20) freeze out and neutrinos decouple. The neutron-to-proton ratio 
departs from equilibrium too, with a value n n /n p = exp (Q/l MeV) ~ 1/6. After 
this time, n n /n p decreases because of the /3-decay of the neutron 



p + e -H- n + v e 
p + + u e -H- n + e + 



(4.20) 
(4.21) 



keep the neutron-to-proton ratio equal to its equilibrium value 




(4.22) 



n — > p + e + v e . 



(4.23) 



At this point, the combination of nucleons experiences a break. All the main BBN 
nuclear reactions involve a deuterium nucleus to be combined with other nuclei. 
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Temperature [MeV] 




Figure 4.6: The evolution of primordial nucleosynthesis, after a PArthENoPE si- 
mulation [148] with r/io = 6.23 and N e g = 3. All the plotted yields but the 
proton and the 4 He ones are normalised to hydrogen, i.e. ni/n p . The ratio n p /nB 
is plotted for protons, while in the helium case we plot the mass fraction Yp = 
4riHe/^B) where tib is the baryonic number density. 



The reaction rate depends on the cross section and on the target abundance, so 
a certain critical deuterium reservoir has to pile up to proceed on. Deuterium is 
produced through the reaction 

p + + n -)> D + 7 . (4.24) 

However, because of 77, the number of photons with energy greater than the binding 
energy of deuterium Bo = 2.22 MeV is easily greater than the number of deuterium 
nuclei, and the inverse reaction is favoured. When the photon temperature has 
decreased below T ~ 0.3 MeV, the abundance of deuterium can finally start to rise 
- as we can see in figure 4.6 — and below T ~ 0.1 MeV the deuterium critical 
density is reached: the deuterium bottleneck opens and heavier nuclei like 3 H and 
helium isotopes can be synthesised. The reaction network is 

D + D^ 3 H + p + (4.25) 
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D + D <->■ 3 He + n 
3 H + p + 4+ 3 He + n 
3 H + D^ 4 He + n 
3 He + D — )> 4 He + p' 



,+ 



(4.26) 
(4.27) 
(4.28) 
(4.29) 



while the production of 7 Be and 7 Li involves the combination of 4 He with 3 He or 
tritium. In the lapse of time till the deuterium bottleneck opening the neutron-to- 
proton ratio decreases to n n /n p ~ 1/7. Because helium is the most tightly bond 
among the light nuclei, almost all the available neutrons are combined into 4 He. 
After the nuclear reactions freeze out around T ~ 30 keV, only very few neutrons 
are left out of 4 He nuclei, mainly in form of deuterium. 

The outcome of standard BBN depends basically on the baryon-to-photon ra- 
tio and on the effective number of neutrino species. The former measures the 
nucleon density, while the latter essentially affects the expansion rate, so these 
two numbers determine the effectiveness of the nuclear reactions rates. Again, 
through a set of Boltzmann equations it is possible to predict the yield of the 
different nuclear species after BBN. Comparing the results of BBN simulations 
with the primordial elemental abundances which are measured thanks to some 
astrophysical observations, it is possible to have an estimate of rj and N cS . The 
overall agreement of the predicted light-element abundances with observations im- 
plies 5.1 < 7710 < 6.5 at 95% C.L. [22], where r] 10 = 10 10 r]. The independent 
determination derived from the CMB temperature fluctuations, r]i = 6.23 ± 0.17, 
agrees rather well with the light-element abundance estimate. A recent estimate 
at 95% C.L. of the number of effective neutrino species from elemental observa- 
tions gives N cS = 3. 68^0^0 and N e s = 3.801q;7q respectively for neutron lifetime 
T n = 885.4 s and r n = 878.5 s [149]. Again, these estimates agree with CMB ones. 
In the following, we will use ^ BBN and r] CMB to refer to the baryon-to-photon ratio 
estimates coming respectively from elemental abundance measurements and from 
CMB analysis. 

The expansion of the universe influences the weak interactions freeze out and 
the deuterium bottleneck opening, and thus the value of the neutron-to-proton 
ratio. Since almost all the neutrons belong to helium nuclei, 4 He was in the past 
the favourite indicator of the presence of extra degrees of freedom in the very early 
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stages of BBN. It is current practice to express the abundance of primordial 4 He 
through its mass fraction 

Y p Hnn/n ) 

n p + n n l + {n n /n p ) 

The value Y p can be estimated from an extrapolation to zero metallicity of the 
measured 4 He content of metal-poor extragalactic HII regions. The systematics 
of the measurements have caused the estimate of Y p to vary significantly over the 
years. We are therefore called for extreme caution when quoting this bound. We 
shall be conservative and adhere to the proposal made in [150]. Here the authors 
set a robust upper bound on Y p based on the assumption that the helium content 
is an increasing function of the metallicity of the cloud. They find 

Y p < 0.2631 at 95% C.L. . (4.31) 

Compared with the standard BBN picture, a higher value of the baryon-to- 
photon ratio implies an earlier opening of the deuterium bottleneck, and interme- 
diate nuclei — like D itself — are processed for a longer time. From the abundances 
of D, 3 H and 3 He we can obtain the value of rj. However, tritium is unstable and 
decays into 3 He, plus an electron and an antineutrino. The estimate of the pri- 
mordial 3 He from astrophysical observation is very difficult, since the only data 
come from the solar system and high metallicity HII regions in our galaxy [22]. 
Moreover, the matter is complicated by the 3 He processing inside stars. Thus the 
abundance of deuterium is the remaining observable. Unfortunately, also the data 
about D are scarce. At the time of the analysis carried on in [1], we only had 
a reliable estimations from 7 high redshift low metallicity clouds absorbing the 
light of background quasars [151], which are plotted in figure 4.7 as a function of 
the H column density. Other two observations recently joined these seven mea- 
surements [152]. The results of these estimations agree well at first glance but 
there is a scatter of the measurements beyond the expectations from the quoted 
systematics. The PDG quotes D/H= (2.82 ± 0.21) x 10~ 5 for the primordial ratio 
of deuterium to hydrogen, where the error has been enlarged to account for the 
still unexplained scatter. Such a large scatter in the measurements of the D abun- 
dance may be a signal of some not yet understood processing of D inside these 
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Figure 4.7: Measurements of primordial deuterium-over-proton ratio D/H. The 
grey bands represent the 95% and 99% range for the observed D abundance. It is 
derived from the 7 high redshift Ly-a clouds shown as a function of the hydrogen 
column density [151]. 



high-redshift clouds. Whether this is the case or not, there are no known astration 
processes that increase the D concentration, so the primordial deuterium yield 
should be larger than this estimation. Any measurement of deuterium constitutes 
therefore a lower limit to the primordial value. In order to be conservative, in this 
study we use 

D/H| p > 2.1 x 1(T 5 . (4.32) 

The dark (light) grey band in figure 4.7 is the experimental 95% (99%) C.L. region. 

Also the 7 Li abundance could be in principle used to determine the baryon- 
to-photon ratio, but it requires a value of r\ rather different from the deuterium 
estimate. However, WMAP measurements of rj agree very well with the deuterium 
value, while the 7 Li one lies 4-5<x out [153, 154]. Because of this primordial lithium 
problem, we are not going to use the 7 Li abundance to constrain axion and ALP 
decays since at the moment observations do not agree with standard BBN predic- 
tions. 

The agreement between independent estimates of rj from BBN calculations and 
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Figure 4.8: D/H as a function of m a for C 7 = 1.9 (Hadronic axions). The width 
of the red band represents the la uncertainty of the CMB determination of r/. 
To the right of the break, axions are treated as being in LTE, to the left they are 
assumed to be decoupled during BBN. The grey bands represent the 95% and 
99% range for the observed D abundance as in figure 4.7. 



WMAP measurements provides one of the most beautiful tests of standard cos- 
mology. A word of caution is however in order since the WMAP value depends 
on cosmological priors such as the spectral index of primordial fluctuations. The 
quoted value stems on a scale-free power-law which we carry as a further assump- 
tion. 



4.2.1 Axions and BBN 

In our decaying pseudoscalar scenario all the BBN parameters can be heavily 
perturbed. In the previous section we have seen the effect on N e g, and here we will 
show how also the concordance between the BBN and CMB values of rj is disturbed. 
The baryon abundance is diluted by the axion decay, as the photons injected in 
the bath increase n 1 suppressing 77. If the decay takes place after BBN but before 
recombination, the concordance between ?7 BBN and 

V CMB ig 

compromised. In our 
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scenario is 



7/ 



BBN 



CMB 



(4.33) 



Taking for granted the WMAP 7 measurement of rj, it means that during the 
BBN time the baryon to photon ratio was enhanced by the factor given by the 
black curves depicted in figure 4.1. This is the main effect influencing BBN. The 
primordial deuterium abundance is the key observable. The increased expansion 
rate plays only a sub-leading role in the axion case. Moreover, here the decay-in- 
equilibrium case provides a dilution factor S2/S1 = 3/2 which is large enough to 
obtain an elemental abundance completely different from the observed one, and 
the equilibrium decay can be excluded too. The bound vanishes for early enough 
decays, which do not perturb too much the standard BBN picture. 

In order to quantify our arguments we have modified the publicly available BBN 
code PArthENoPE [148] to include the effects of axions, taking into account their 
impact on the Friedmann equation and the modified densities of different radiation 
species. Our results are the main outcome of [1] and are shown in figure 4.8, where 
we plot the D/H as a function of the axion mass. For m a > 20 keV we treat 
axions as being in thermal equilibrium throughout BBN. For m a < 10 keV we use 
their abundance from our numerical freeze-out calculation, assuming that they are 
decoupled during BBN. We do not treat the intermediate case, leaving a gap in 
the predicted deuterium yield as a function of m a that is shown as a red band in 
figure 4.8. To calculate this curve we have adjusted the baryon abundance such 
that in the very end it matches the CMB-implied value. Its la range is represented 
by the width of the red band. 

We compare the predicted D yield with the measured value, and the results 
are plotted in figure 4.8. As we see in figure 4.8, the presence of axions reduces 
D/H, so our bounds appear to be conservative regarding the possibility of some 
unknown stellar production process. 

The deuterium abundance is reduced below its 2a observation if m a is below 
300 keV. Therefore, BBN constrains axions to have masses 



Axions with mass above this limit have almost completely disappeared from the 
thermal bath before they can affect BBN. In this case the predictions approach 



m a > 300 keV . 



(4.34) 
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standard BBN. This mass upper bound closes the red Hot DM band in figure 1.2. 
Since this bound corresponds to axions that are always in thermal equilibrium, even 
only via the Primakoff and inverse decay processes, it also applies to non-hadronic 
axions which would interact more strongly. 



4.2.2 ALPs and BBN 

Again, the phenomenology of ALP decay leads to a broader sample of situations 
and requires a scan of the parameter space. In order to numerically evaluate the 
impact of decaying ALPs in the BBN predictions, we have used a BBN code that 
includes the modified cosmology driven by ALP decays computed with the tools 
of the previous sections. We have written a simple BBN code in Mathematica to 
compute the primordial abundances of D, 3 He, 4 He, 7 Li and 7 Be. We have used 
the minimal reaction network relevant for rj ~ )j CMB and N c g ~ 3 as detailed 
for instance in [155-157]. When facing standard ALP-less cosmology, our results 
are in very good agreement with standard BBN calculations obtained with the 
KAWANO [158] or PArthENoPE [157] codes, given the theoretical and experimental 
uncertainties, which gives us confidence in our results. This allows us to easily 
compute the outcome of BBN when ALPs have a non trivial role during BBN 
itself, and gives the right trend when ALPs decay much later than BBN, enhancing 
enormously the value of ^ BBN with respect to ^ CMB . In the latter case, the minimal 
reaction network we used might be not sufficient to have a precise BBN outcome, 
but anyway it gives a first order approximation that leads clearly to exclude a 
too large ^ BBN . The effect of relic particle decays on the abundances of light 
elements predicted by BBN is largely discussed in the literature [159-161], usually 
considering very large mass particles, with m <C GeV. In our analysis instead we 
focused also on lighter masses, in the keV and MeV range. 

The impact of ALPs on BBN, which was discussed in [2], depends strongly on 
the m^, in particular whether ALPs are heavy enough to induce electromagnetic 
or hadronic cascades. In the following these cases are discussed separately. 
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Small masses 

If the ALP mass is smaller than a few MeV, the decay products cannot induce 
nuclear reactions and their effect on BBN is only indirect. The injected photons — 
and perhaps a small amount of electron-positron pairs — dilute both the neutrino 
and baryon densities. The impact then depends on whether the decays happen 
before or after BBN. 

For decays happening after BBN, the injected photons heat the bath, decreasing 
the baryon to photon ratio rj. The outcome of this high-?7 BBN scenario is the same 
of the axion case: intermediate nuclei like D or 3 He are more easily consumed, and 
the final abundance of heavier nuclei like Li increases. 

However, the effect of ALP decay can have a strong effect on the 4 He yield 
too. In the late ALP decay scenario, the bottleneck opens earlier because of the 
enhanced 77, so neutrons have less time to decay, enhancing the final 4 He yield. 
But this is not the only effect. The presence of ALPs makes the universe expand 
faster, which has two additional implications. It induces an earlier freeze out of 
thep -H- n conversion reactions, and therefore a larger n abundance. Moreover, the 
time between this freeze-out and BBN is even more shortened, and therefore the 
amount of neutrons that decay is lower. Thus three mechanisms are responsible 
for the enhancement of the 4 He yield. 

We have r/ BBN = ^ CMB if ALPs decay before BBN, i.e. before the opening of the 
deuterium bottleneck. The main trends mentioned before disappear. However, the 
ALP decay can still modify BBN indirectly, if it happens between the freeze-out 
of weak interactions and BBN. There are three effects that we should take into 
account. First, ALPs are present during the freeze-out of the p -H- n reactions, 
so the n abundance is in principle larger because of the faster expansion. Second, 
when ALPs decay they reduce N eS , as shown in the previous section. The cosmic 
expansion is slow and neutrons have more time to decay. The neutron concentra- 
tion is thus affected in two opposite ways which tends to compensate each other. 
The time that neutrons have to decay depends on how close to BBN the decays 
happen. In fact we find that 4 He grows as the ALP decay happens closer to BBN. 
Only in a small 4 He region at ~ MeV and r ~ 30 s, neutron decay plays a 
role inducing low 4 He. Because of the two opposing effects, the 4 He abundance is 
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therefore not a sensitive indicator for ALPs in this region. Anyway, this scenario 
also implies low D and 3 He and high Li, even if for a completely different reason. 
The slower expansion gives more time for the consumption of intermediate nuclei, 
and the final D yield is again lower than the standard BBN one. 

The predictions for the primordial mass fraction of 4 He, Y p , and the deuterium- 
over-proton ratio D/H in the ALP-decay scenario are shown in figure 4.9. The 
excluded regions are above the thick isocontour. The isocontours very much re- 
semble those of N eS because the outcome of BBN is mostly sensitive to the value 
of r] BBN and therefore to the baryon dilution, which qualitatively follows the same 
logic as the neutrino dilution. ALPs with small mass and fast decay disappear 
from the bath in local thermal equilibrium with photons at temperatures around 
m$. Therefore here BBN only depends on the ALP mass, not on the lifetime, and 
the outcome is the same of equilibrium-decay of axions that provide the bound 
(4.34). When the lifetime is longer, the isocontours are parallel to the lines of 
constant entropy production as in the N cS case. 

The exclusion bound from 4 He is depicted as a purple region in the chapter 
summarising picture 4.10, while the deuterium bound is shown in red. In the light 
masses range is deuterium that provides the main exclusion bound. The bound 
corresponds to ALPs that inject a a fraction of order 10% of the total entropy 
in the electromagnetic bath after BBN. Under this circumstances the baryon and 
neutrino dilution is 0(1), and, since deuterium is the most sensitive observable 
to r] BBN , it is also the most constraining argument. Note that the 4 He abundance 
depends only logarithmically on ^ BBN while D/H oc 1/(77 BBN )~ L6 [153]. 

Large masses 

The decay of large mass ALPs affects BBN through N eS and rj, which continues 
along the trend of the analysis for small masses. The additional ingredient in 
this case is the direct dissociation of nuclei due to energetic enough decay pro- 
ducts, which can strongly modify the predictions of BBN. We have depicted the 
areas where these effects are relevant as grey regions in figures 4.9. Moreover, 
the creation of pions in hadronic cascades has also a peculiar way of affecting the 
primordial nucleosynthesis. 
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Figure 4.9: Isocontours of the primordial abundance of deuterium normalized 
to protons D/H (left) and helium Y p (right) in the decaying ALP cosmology, as 
a function of the ALP mass and lifetime. 



We first comment on photo-dissociation, since the main ALP decay channel 
we are analysing is the two-photon one. The photo-dissociation of deuterium - 
which requires a threshold energy E t h = Bo = 2.22 MeV - - of course reduces 
its abundance further, making our bounds even stronger. On the other hand, 
the photo-dissociation of 4 He can create deuterium, which reverses the trend of 
our constraints. This requires higher-energy photons since the threshold energy is 
-Eth = 19.81 MeV. We have followed the methods of [162] to check in which region 
of parameter space this phenomena can change our predictions. At high densities 
and temperatures, high-energy photons interact very fast with the thermal bath 
creating electromagnetic cascades. The energy injected is rapidly redistributed 
and the resulting cascade spectrum is [163] 

if E < E x 

n T) 

1 (E) = tK (^Y if E x < E < E c , (4.35) 



dE 




otherwhise 
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where the normalisation constant Kq is determined equating the total and the 
injected energies. The spectrum features a knee and a very sharp cut-off at high 
energies which depend on the plasma features and temperatures, E x ~ ml/(80T) 
and E c ~ m 2 e /(22T). At high T, this cut off lies below the photo-dissociation 
threshold of nuclei and these effects are negligible. Therefore, requiring E c (T d ) > 
E th and converting the decay temperature in seconds, we see that ALP decay 
cannot dissociate deuterium for r < 3000 s and Helium for r < 3 x fO 5 s. A 
second key quantity is the injected energy per photon in the bath 

Q = m^— = — - . (4.36) 

™7 Z 9*S{- L foj 

The fractional change for a nuclear abundance is [162] 

5Xi ( X T 



where (E) is the average energy in the electromagnetic cascade and fr^p is the 
relative strength of the cross section for the photo-destruction of the target T into 
the product P compared with the Thomson cross section. Taking the deuterium 
case as an example, a fractional change of 1/2 is provided by [162] 

.prod o 9yln -ll rp yf 5000 W 5.0xl0- 4 \ /^qX / E th \ 1/2 / T X ^ 

(4.38) 

if production dominates, while if destruction does it is [162] 

cr-, xl o- 8Ge v(^) (T )(^) 1/2 (^)-. (4 , 9 ) 

In our analysis, the ALP thermal origin makes the fractional change for D and 
4 He always too strong, and the regions where photo-dissociation is allowed can 
be safely excluded unless some extra dilution of the ALP population is taken into 
account. 

For ALP masses above a few GeV, ALP decays produce radiatively quark- 
antiquark pairs that will hadronise. Hadronic cascades can dissociate nuclei in a 
similar fashion of the electromagnetic ones, if they are happening after a typical 
time of r ~ 10 2 s [164-166]. For earlier decays, the electromagnetic interactions in 
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the plasma make the hadrons lose rapidly their energy and the nuclear destruction 
is suppressed. The effects of hadronic cascades are again necessarily dramatic 
because of the large ALP relic abundance. We consider extremely unlikely that the 
combined effect of nonstandard BBN with the post-BBN processing gives similar 
results to standard BBN. So we exclude all the regions where electromagnetic and 
hadronic cascades play a role, see figure 4.9. In the literature, cascade constraints 
are usually presented in terms of plotted in function of lifetime r. In the case in 
exam, it is possible to refer to this representation almost directly from figure 4.9, 
considering that the ratio n^/n^, given by equation (2.54), is constant in all the 
plotted parameter space in which cascades play a role and does not change much 
outside. 

A crucial difference between electromagnetic and hadronic cascades lies in the 
peculiar effect of pions produced in the latter. If is larger than twice the 
charged pion mass = 139.57 MeV, the decay channel — > 77r + 7r~ opens up. 
Even if its branching ratio is very small, the abundance of relic ALPs is such that 
a huge amount of pions — compared with that of the present nuclei — can be 
produced. If the decay happens before or during BBN — typically r < 100 s - 
the universe is dense enough to make pions induce neutron-proton interconversions, 
7r + +n — > n°+p + and ir~+p + — > ir°+n, before decaying [167]. The second reaction 
is favoured because of the Sommerfeld enhancement and the typical overabundance 
of protons over neutrons at T < Q. These reactions will therefore tend to increase 
the neutron to proton ratio. They can do it much even more drastically than 
the mere presence of the ALP during the p ■<->■ n freeze out commented in the 
previous subsection. The higher neutron abundance would aid the heavy element 
production but it also increases D/H, which is the most important effect. Since 
almost all neutrons end up in 4 He nuclei taking protons with them, a higher initial 
neutron abundance yields a smaller final proton abundance and thus a larger D/H 
ratio. If n n /n p ~ 1 at the onset of BBN, all protons end up in 4 He and D/H would 
be arbitrarily large! 2 

We have included the effects of pions in our BBN code following reference [165]. 



2 Of course including neutron decay during BBN would still give a finite, albeit very large, 
result. 
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In figure 4.9 we can see that the low D/H trend of low mass ALPs is drastically 
changed when the mass gets above 2m n + and the effect on 4 He gets strongly 
boosted when crossing this boundary. The effects of pions are hampered if the 
ALPs decay very early for r < 1CT 2 s. In this case the electroweak reactions 
p + + e~ -H- n + u e can still re-establish the n n /n p equilibrium. Also very late decay 
t > 10 2 s, does not affect BBN through pions injection, because ns fail to interact 
before decaying. In this last tiny region our results cannot be taken quantitatively 
on trust, since we have not taken into account the possibly ineffective slowing down 
of pions after e + e~ annihilation. 

The region of the parameter space where cascades are protagonists are depicted 
in pink in figure 4.10. Again, the purple region is excluded by 4 He overproduction. 
Together with the neutrino dilution limit this is the best bound for the > 2m n 
region of ALP parameter space. The parametric dependence of the bound is 
the same of the iV eff limit. Like in the axion case, short-lived ALPs that decay 
before BBN can not be constrained. In this large mass region, it corresponds to 
r < 10" 1 - 5 s. 



4.3 Effect of other couplings 

A few words are in order to consider the effects of other couplings, besides the 
two-photon one, in the ALP case. Additional couplings imply in general a higher 
interaction rate and therefore a later decoupling and a larger number density. 
Of course, the more abundant is a population the more important are the conse- 
quences of its decay. However, in our bounds these considerations are only partially 
true. Since we are not directly taking into account possible degrees of freedom be- 
yond the standard model, a population of ALPs with additional couplings would 
be more abundant only if the new decoupling temperature is lower than E EW , 
otherwise the dilution factor does not change, because g*(T > E EW ) = 106.5. 

For a generic PNGB a CP-odd coupling to GG and derivative couplings to 
fermions are allowed, 

£*ff> = Wffl^r , (4-40) 

the latter even flavour nondiagonal. The coupling to fermions (4.40) allows the 
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Figure 4.10: Bounds on early ALP decays from deuterium underproduction (D, 
red), helium overproduction (He, purple), helium photodissociation (7He, pink), 
hadronic cascades (Hadr, pink) and neutrino dilution (N e s, yellow). 



decay (f> —iff at a rate 

3/2 



CffY i. m f + m rf m <t> L ( m f + m f^ z 



rrif — rrifi 

m<t> 



16tt v \ 

(4.41) 

which is suppressed with respect to the two-photon decay for small fermion masses. 
Writing = C 7 a/(27r/^) the <p — > ff can dominate only in an interval near the 
kinematic threshold 1 > (rrif +m'f)/m ( f ) > aCy/AirCffi. For ALP masses above 
few GeV, the coupling to gluons allows the ALP decay into two gluons at a rate 

r</>-> 99 = 8(C 99 /C 7 ) 2 r 9 i_ >77 . 

If the ALPs are cosmologically stable, the bounds from direct detection of ALP 
decay photons and the DM overproduction — which are described in the previous 
chapters — still hold. These limits depends mainly on the ALP lifetime and abun- 
dance. The pre-decay abundance can only be increased by additional couplings, 
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but the area constrained by these limits lies all in the Tf G > -Eew region, where 
is constant according to our assumption on the primordial particle content. For a 
given photon coupling, the lifetime is shortened if more decay channels open. Thus 
the limits on the long lifetime range are only mildly affected by the abundance in 
the m^-T parameter space, while in the m^-g^ plane they have to be additionally 
adjusted according to the different relation between lifetime and photon coupling. 

The short lifetime region corresponds to the constrains described in this chap- 
ter. At low masses, the deuterium and helium bounds come from ALPs in thermal 
equilibrium with the bath. Clearly, adding more couplings between the ALP and 
SM particles we cannot avoid these bounds. In the intermediate mass region 300 
keV < < 2 m M , where m M = 105.7 MeV is the muon mass, the D/H and N eS 
bounds follow from the dilution of baryons and neutrinos with respect to photons. 
These bounds apply to ALPs decaying into photons or electrons, and indeed we 
have not made a difference between the two in our equations. Again, in the m^-r 
plane these bounds suffer only a mild change due to the different abundance. When 
we translate the bounds in the m^-g^ plane they will show worst if the decay into 
two electrons dominates than if we only consider the two photon coupling, because 
what is important for the bound is the total decay rate. The lower bound on 
reduces by a factor ~ (2nC ee m e ) / (aC^m^). 

Since the direct decay into neutrinos is suppressed by ~ [m v j 'm^) 2 , an ama- 
zingly tiny number, the bounds are perfectly valid provided one interprets r as 
the total lifetime, not only due to the two photon decay channel. Of course, this 
is valid unless one considers sterile neutrinos with m u ~ m^. Then in this case 
neutrinos have a strong tendency to constitute too much DM. A way to avoid 
this is to make them decay into a SM neutrino plus a photon, but this produces 
entropy and we expect a similar, slightly smaller, bound from D/H in this case. In 
these models the low N cS tendency is reversed since the sterile neutrinos produce 
neutrinos in its decay. 

If 2m n + > > 2m M we have a somewhat different scenario where the ALP 
tends to favour the — > fi + fi~ decay. The upper limit on r in this region comes 
from having too low N eS already before BBN. But if the decay into muons domi- 
nates we will rather have a high N cS . In this case the amounts of energy released in 
electrons and in neutrinos by muon decay ji — > euu are similar. Since data favours 
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values larger than the standard N e g = 3, the N e e bound will relax somehow. Any- 
way, we do not expect them to disappear, because ALPs can still produce too 
many neutrinos. Also in this case the bound on deuterium should come from a 
too high D/H, which is less conservative a constraint. In any case the bound from 
He will stay since it mainly comes from a high ry BBN and the ALP contribution to 
the expansion at the freeze out of p ■<-)■ n weak reactions. 

Finally, for > 2m w + the most stringent bound comes from 4 He overproduc- 
tion due to the presence of charged pions before BBN, enhancing the neutron to 
proton ratio. As we commented, this bound does depend very little on the details 
and branching ratios of the ALP since only a minimal number of pions are sufficient 
for a drastic change. Therefore we expect it not to change very much. However, 
when quoting this constraint in the m^-g^ plane this bound would display in a 
lower position than in the case where only the two-photon coupling is considered. 
Only in this region the coupling to two gluons can affect the ALP decay and will 
certainly increase the pion multiplicity of the decay making the bound on r slightly 
better. The decay into muons can dominate if m<j, is not too far from 2m M and all 
said in the above paragraph holds. It appears that the helium bound will still be 
the most relevant in this case. 



CHAPTER 5 



Summary and conclusion 



The axion is a side product of the elegant solution of the strong CP-problem 
proposed by Roberto Peccei and Helen Quinn in 1977. The different realisations 
of the Peccei-Quinn idea produced a variety of axion models, which are tested 
in particle physics laboratories, and challenged in astrophysical and cosmological 
observations. Moreover, it seems that particles with similar characteristics to the 
axion could arise in several extensions of the standard model of particle physics. 
In particular, string theory seems a fertile environment that can provide plenty 
of these axion-like particles. One day, the so craved experimental test of string 
theory could finally come from ALP-related observations. 

In this dissertation we have depicted the constraints that cosmology puts on 
the existence of such particles. In chapter 1, we started with a presentation of the 
axion theory, together with a brief motivation for ALPs. Then, we introduced the 
general limits on the parameter space of axions and ALPs. We listed the direct 
experimental tests and the astrophysical observations that exclude the existence 
of these pseudoscalar for determined choices of mass and couplings. But to give 
new hope, we also reviewed some astrophysical problems that could be solved by 
pseudoscalar particles in two different ranges of the parameter space. 

Motivated by the possible solution of particle and astroparticle problems, we 
continued our analysis in chapter 2, where we treated the possibility for a primeval 
population of pseudoscalars to arise in the early universe. This enquiry sets the ba- 
sis for the subsequent discussion about cosmological limits, and most importantly 
provides the motivation to consider axions and ALPs as a constituent of the dark 
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matter of the universe. So the existence of these particles can provide the solution 
of a further problem of modern physics. In this chapter we illustrated a very rich 
phenomenology. Axions and ALPs can be created by thermal interactions with 
the particles of the primordial plasma. But the non-thermal creation — via the 
so called realignment mechanism — is even more interesting, since it involves the 
physics of spontaneous symmetry breaking and phase transitions applied to the 
primeval universe, and may remarkably provide the explanation to the dark mat- 
ter. A first limit on the parameter space can be put at this stage, if too much dark 
matter is produced. 

Also the cosmo logical stability of pseudoscalars is discussed in chapter 2. Their 
peculiar two-photon coupling provides them with a decay channel. Therefore, a 
pseudoscalar population can decay in photons or be reabsorbed by the primordial 
thermal bath, if this electromagnetic interaction is active. We briefly treated also 
the role primordial magnetic fields could have in this process, to further develop 
the discussion in the appendix A. 

The decay of a population of particles during the early epochs of the universe 
can have dramatic consequences and leave an indelible imprint on cosmological 
observables. This is the topic of the second part of this study. We have divided 
the limits into two sets, according to the epoch of the decay. Both of them are 
dealing with the influence of photons produced by the pseudoscalar decay on the 
successive evolution of the cosmos. 

In chapter 3 we collected the limits related to late decay and to photon de- 
tection. These could either be the CMB photons, which could have imprinted in 
their spectrum the distortions caused by the decay products, or the decay photons 
themselves, which could have freely travelled towards our telescopes. Observations 
of photon spectra leave very little room to pseudoscalar decay — especially in the 
CMB case — and the decay of a whole population can be safely excluded, unless 
it happened when electron-photon interactions were active and able to thermalise 
rapidly the injected photons. In this chapter we constrained rather long lifetimes, 
since the distortion of the CMB requires the decay to happen just before the re- 
combination era, which occurred when the universe was roughly 400,000 years old. 
The direct observation of the decay products excludes cosmologically stable ALPs 
and axions up to r ~ 10 24 -10 28 s. A section in this chapter is dedicated to the 
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ultraviolet radiation that can be eventually emitted in late decays. The universe 
is very opaque to ultraviolet light. This kind of photons can not appear in direct 
observations, but they can nevertheless be constrained because of their effect on 
the ionisation history of the universe. 

Early decays that do not affect the CMB spectrum can be constrained by the 
arguments of chapter 4. The entropy transfer and increase due to the decay of a 
relic population have the leading role in the discussion. We first discussed the case 
of a population that decays in local thermal equilibrium. Entropy is conserved and 
just transferred to the species in thermal contact with the disappearing popula- 
tion. The effects of an out-of-equilibrium decay can be more dramatic, especially 
if the pseudoscalar population is dominating the energy density of the universe 
before the decay. A large amount of entropy is created and transferred to the 
photon bath. In both cases, the subsequent evolution towards thermal equilibrium 
makes the temperature of the species in thermal contact with photons to increase 
relatively to the decoupled ones. In this sense we defined this event as a dilution of 
the decoupled species, the best example being neutrinos. Solving numerically the 
set of Boltzmann equations that describes the evolution of pseudoscalar, photon, 
electron and neutrino populations, we calculated the effect of the decay on neutri- 
nos. We then compared our neutrino dilution scenario with cosmological neutrino 
observations. Through CMB multipole analysis and LSS survey it is possible to 
measure the number of effective neutrino species N cS . Present data prefer iV eff > 3, 
which works against the cosmological dilution of neutrinos. After this observation, 
we were able to put some limits on pseudoscalar parameter space. 

The outcome of primordial nucleosynthesis is also influenced by the pseu- 
doscalar decay. Both directly, if pseudoscalars are massive enough to inject decay 
products energetic enough to break nuclear bonds, and indirectly, through the in- 
fluence on the number densities of baryons and neutrinos relative to the photon 
one. This topic is also part of chapter 4. The primordial yield of deuterium is very 
sensitive to the baryon-to-photon ratio, which is heavily perturbed in our scenario 
of early pseudoscalar decay. It indeed provides the most restrictive cosmological 
lower bound on the axion mass, m a > 0.3 MeV. Also ALPs are severely constrained 
by BBN, since they affect the baryon-to-photon ratio, N cS and can generate de- 
structive electromagnetic and hadronic cascades. The production of pion cascades 
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Figure 5.1: Cosmological ALP bounds in the m^-g parameter space. The la- 
belling is described in the text. 

in the decay of > 2m n ALPs has a peculiar effect on the primordial 4 He out- 
come, which we constrained too. Chapter 4 concludes with some considerations on 
the changes in the cosmological bounds that further ALP couplings could provide. 
The conclusion here is that cosmological limits, although slightly modified, are 
very solid. 

In our summary plots 5.1 and 5.2 we collected all the bounds arising from co- 
smology, considering a thermal origin for the primordial pseudoscalar population. 
While figure 5.1 is meant to be the reference for the particle physicist, since the 
ALP parameter space is plotted in function of the ALP mass and coupling 
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#0, figure 5.2, where the parameters are and the lifetime r, provides a better 
physical understanding. In this picture we can better follow the relation of the 
limits with the sequence of events in the history of the universe. The cosmolo- 
gical bounds critically depend on the ALP lifetime while they have only a milder 
dependence on other parameters. In figure 5.1 they have indeed the characteristic 
slope oc m^ -3 / 2 . As a reference, we plotted the line along which ALPs have the 
lifetime equal to 10 17 s, i.e. the age of the universe. To warn the reader about our 
assumptions, we also plotted the lines for freezing out temperature Tf Q equal to 
the electroweakscale -Eew an d the Planck scale m-p\. Most likely, Tf Q > mpi makes 
no physical sense, and for Tf Q > -Eew the actual abundance could be lower than 
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that used in our calculations, since we have not taken into account other particles 
besides the ALP and the SM ones. 

Ordered by decreasing lifetimes, the excluded regions are: 

• DM — if ALPs are cosmologically stable and behave as dark matter they 
should not exceed the DM fraction measured by WMAP. 

• Optical, X-Rays, 7-Rays — photons produced in pseudoscalar decays 
inside galaxies would show up as a peak in galactic spectra that must not 
exceed the known backgrounds. 

• EBL — photons produced in pseudoscalar decays when the universe is tran- 
sparent must not exceed the extragalactic background light. 

• Xi on — the ionization of primordial hydrogen caused by the decay photons 
must not contribute significantly to the optical depth after recombination. 

• CMB y, n — if the decay happens when the universe is opaque, the de- 
cay photons must not cause spectral distortions in the CMB spectrum that 
cannot be fully rethermalised. 

• EM, Hadr showers — the decay of high mass ALPs produces electromag- 
netic and hadronic showers that must not spoil the agreement of big-bang 
nucleosynthesis with observations of primordial nuclei. 

• 4 He, D — the ALP and axion decays produce photons that dilute the baryon 
and neutrino densities, whose values affect the outcome of BBN, in particular 
the deuterium and 4 He yields. Again, this dilution should not compromise 
BBN. 

• N e flf — the neutrino density must not disagree with the value measured by 
WMAP and other large-scale-structure probes. Currently, data points to an 
effective number of neutrinos iV eff greater than 3, which is disfavoured in the 
cosmology of decaying pseudoscalars. 

We have seen how cosmological observations can exclude a large part of the 
pseudoscalar parameter space. These limits are solid, and the prospects are good, 
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since the amount and quality of cosmological observations is steeply rising. This is 
very important, since the constrained region lies in a part of the parameter space 
presently inaccessible to direct experimental tests. 

The investigation on the effects of axions and ALPs on astrophysical and co- 
smological observables must proceed further, because the detection of one of several 
striking signatures could lead to their discovery. The smoking gun could be hidden 
in the white-dwarf evolution or in the structure of galaxies, maybe influenced by the 
dark matter in a Bose-Einstein condensate, or in some different phenomena that 
we still have to analyse. But after the eventual and desired discovery, the direct 
detection and laboratory experiments will have the task of measuring precisely 
the new particle characteristics. In the near future, direct detection experiments 
like the haloscopes and the helioscopes will finally reach the sensitivity to test the 
axion hypothesis and to explore two phenomenologically important regions of the 
parameter space, and maybe some long standing questions will find an answer. 

Understanding the ultimate theory behind the laws of Nature is the final pur- 
pose of physics. The large hadron collider is pushing the knowledge on particle 
physics towards unexplored energy scales. The recent discovery of the Higgs boson 
gives a bit more confidence about the axion theory, which requires the existence 
of scalar fields and the spontaneous breaking of symmetries at high energy scales. 
But increasing the energy tested by colliders is not the only way to probe the 
physics beyond the SM. If axions and ALPs will be finally discovered, we will 
have the tempting possibility of exploring the physics related to very high energy 
scales through its low energy regime. Still largely unexplored, the low energy 
frontier could hide important novelties and bring to fundamental advances in the 
understanding of Nature. 



APPENDIX A 



Axion-photon mixing 



In the following we follow the treatment of [168] and [169], where the hidden 
photon resonance case was studied. The Lagrangian of the pseudoscalar-photon 

— * 

system in presence of a strong magnetic field B cxt , whose modulus is B ext , is 

C = —F^F^ + + U^d^ - X -m^ 2 - g^B cxt ■ d A , (A.l) 
where we expressed FF/A as the scalar product of the external magnetic field and 

— * 

of the electric field component — d$A of the field strength. From the Lagrangian 
we obtain the equations of motions for the pseudoscalar field and the component 
of the vector potential parallel to the magnetic field A// 

HA,/ - ad A// + g^B cxt d <f) = , (A.2a) 
□0 + + g+B^doA, = , (A.2b) 

while the perpendicular component A± is not affected by the interaction with 0. 
In equation (A. 2a) we have used Ohm's law for the current density, j M = — adoA^, 
assuming a linear response of the medium whose conductivity is o [125]. In Fourier 
space the equations (A. 2) become 

(— uj 2 + k 2 — iuia) A a + ig<f,B cxt u(fi = , (A. 3a) 

(-co 2 + k 2 + m/) + ig^B cxt uA// = . (A.3b) 

Assuming the quantity x — \/ 94>B cxX ,oj to be negligible with respect to the other 
energy scales of the problem, this system of equations can be diagonalised by 
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means of the redefinitions 

A,, = A/, - (A 4a) 

iujo — m<f, 2 

- mB^ k (A4b) 

tcoa — 

The new states A// and are decoupled, and if the mixing parameter 

2 



Xeff 



icua — m^ 2 



(A.5) 



is very small, they are almost photon and pseudoscalar states. They can oscillate 
one into another, if ~ u. Their dispersion relations are respectively 



u 2 - k 2 A// = -iua + 0( X 2 ) , (A.6a) 

.2 l2 _ _2 , <n(,.Z 



u, 2 - A;J = ml + 0( X 2 ) , (A.6b) 

where the 0(x 2 ) factors are complex. In the case of zero mode DM scalars, = 0. 
The original states A// and can be obtained inserting the A// and solutions 
into the (A.4). 

For the conductivity of the medium we use the relation [3] 



«cru; = cjp g 1" Th = ^ 7 (-t ) + iojD(u, 1 ) , (A.7J 

1 + (wr col i) 



that provides an interpolation of the classical approximation and the quantum 
regime. In this formula, cup is the plasma frequency and Txh = o"Th^e is the 
photon absorption width due to Thomson scattering off electrons, while r co u is the 
average time between electron collisions, which sets the time scale. If we call p 
and E the momentum and the energy of electrons, the plasma frequency is [46] 



2 4a [°° p 
cup = — / ap- 

n Jo 



exp (E/T) + 1 
We next define the squared plasmon mass 



E 3 \EJ 



(A.8) 



. <iiian e /m e if T <C m e 
m 2 = Re (iau) = { ' . (A.9) 

\caxT 2 if T > m e 
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The damping factor 

D(u, T) — Im (ia) (A. 10) 

expresses the net-rate at which photons are created or absorbed. The full ex- 
pressions for these quantities can be found in [168]. Using these definitions, the 
oscillation probability (A. 5) becomes 



xlfi 



ig<f,B ext uj 



icoa — m^ 2 



(g^B ext u) 2 



(m 2 - ml) + (ooDY 



(A.11) 



Because of this magnetic field-mediated mixing with the photon, the pseudoscalar 
field has the same interactions of the photon, but with a rate suppressed by xls ■ 111 
particular, it can be absorbed or emitted during Thomson-like interactions 0+e -H- 
7 + e and annihilations- like events + 7 -h- e + + e~. If thermal contact is regained 
by the pseudoscalar condensate, its distribution tends towards equilibrium. The 
amount of pseudoscalars that survives the magnetic field induced evaporation is 
provided by 

n^to) = n^U) exp (- J ° X 2 cS D dtj , (A.12) 

where the subscripts i means quantities measured just after the realignment me- 
chanism ended and are today quantities. We used in the integral the damping 
factor D, that provides the net-rate of absorbed photons. Because we deal with 
temperature-depending quantities, for future convenience we can write the integral 
in the exponential of equation (A.12) as 

f t0 9 „ , f To (gsB cxt uj) 2 DdT , k . 

-/ xlvDdt- / w 9 ; . A.13 

The primordial magnetic field evolves with time, since it is stretched by the 
expansion of the universe. The conservation of magnetic flux implies B cxt = 
B cx t,*{R*/ R) 2 ~ B cxt ^(T*/T B ) 2 , where the subscript * refers to quantities mea- 
sured at the reference time i*, after the creation of the primordial magnetic field 
at T B . 

The temperature dependence of the conductivity (A. 7) is mainly affected by 
T co\\ ~ (^Vy^e) -1 oc T^ 1 . For T ^> uj the conductivity is dominated by the damping 
factor, thus xls ~ (fiV>-B cx t) 2 / D 2 if also ra^/u <C D. The damping factor is 
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D = T/(3na) in a relativistic plasma of electrons [125], and equation (A. 13) 
becomes 

To (g<t>B ext ) 2 D dT g\Bl f T o T 4 mpi dT 



[ l0 ( g4> B cxt y PdT gffi r 
J Tb ' D 2 HT ~ T 4 i T 



T 4 7 Ts 77(3™) 1.66v^(T)T2 T 

■)•(; 



B \ 2 / 106 N 1/2 



10- 10 GeV-V VnG; \g*(T B ) J 10 9 GeV 



(A.14) 



For axion cold DM this is the approximation to use, because the axion mass is 
much lower than the damping factor [125]. Axion cold DM is not affected by 
primordial magnetic field: its two photon coupling strength is much lower than 
10~ 10 GeV -1 , while the present value of the magnetic field B can not be stronger 
than few nG [126,127]. Even if Tb is very large, it cannot affect the condensate at 
temperatures larger than f a , because in this case there is no condensate. Moreover, 
in the case of the hadronic axion, there is no two-photon coupling above Aqcd- 

When the damping factor does not suppress XeS and it is much smaller in 
magnitude than the plasmon mass, we can distinguish three different regimes of 
this quantity, according to the relative magnitude between m 7 (T) and m<f,. The 
resonant regime for m 7 (T) ~ is the most important, since the mixing gets 
enhanced by the cancellation of the mass terms in the denominator. 

At temperatures higher than the resonant one, equation (A. 13) is 

(g^gext^) D I ^fl WJ^[ I > m e 



(rrij 2 



I a 2 (n e /m e ) 2 T' A /m m e 

(A.15) 

which in both cases is dominated by lower temperatures, so towards the reso- 
nance. For the low temperature we used the conductivity provided by a hydrogen 
plasma [125]. The low temperature regime gives instead 

2 f {g*B^)\T/T.f T/a 

(g^B^u) D I ^ r3/ mpi 1 » m e 



m^f + {uDf HT | (g^ext,^) (T/nf T 3/2 me -l/2 a -l T 



(A.16) 

which is dominated by the higher temperatures again towards the resonance. The 
behaviour of the integrand can even better understood from figure A.l, where this 
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quantity is plotted for = 0.1, 1 and 10 eV. Since this function is so peaked 
around the resonance, we can approximate the integral evaluating all the temper- 
ature dependent quantities around T res , expanding ~ m| + drn^/dT (T — T res ), 
and neglecting the rest 



L 



T ° (gjB^u) 2 D dT ^ 
Ti (mi* - m^) 2 + {uD f H T 

T ° (^^cxt.rcs^) 2 D dT_ 

Ti (dm*/dT) 2 (T - T rcs ) 2 + (ojD) 2 H rcs T rcs 

71 (9cf>B ex t,ves) 2 1 



2 H rcs m^ (d\nm? ( /d\nT) res ' 
To obtain the last line we used the approximation 

dx 7r 1 + arctan(y/,2) 



I 



(A.17) 



(A.18) 
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and oj ~ m^. The pseudoscalar population is affected by the resonant magnetic 
field evaporation if the value of the integral (A. 17) is larger than one. This is 
plotted in figure 2.10. 
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